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Abstract 

We consider a nonlinear area preserving Anosov map M on the torus phase space, 
which is the simplest example of a fully chaotic dynamics. We are interested in the 
quantum dynamics for long time, generated by the unitary quantum propagator M. 
The usual semi-classical Trace formula expresses Tr (m*\ for finite time t, in the limit 
H — > 0, in terms of periodic orbits of M of period t. Recent work reach time t <C £e/6 
where tE = log(l/ft) /A is the Ehrenfest time, and A is the Lyapounov coefficient. 
Using a semi-classical normal form description of the dynamics uniformly over phase 
space, we show how to extend the trace formula for longer time of the form t ~ C±e 
where C is any constant, with an error O (ft? ). 
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1 Introduction 

Semi-classical analysis is a fruitful approach to understand wave equations in the regime 
where the wave length is small in comparison with the size of the domain, or with the size 
of the typical variation of the potential, where the wave evolves. In that regime, one shows 
that the evolution of a wave can be described in terms of Hamiltonian classical dynamics 
in the same domain (or with the same potential). For example, the Van-Vleck formula 
(1928) expresses the evolved wave as a sum of the initial wave transported along several 
classical trajectories. Because the wave formalism enters in many area of physics (acoustic 
waves, seismic waves, electromagnetic waves, quantum waves ...), semi-classical analysis is 
an important mathematical tool to understand physical phenomena. Wave equations, and 
more generally Partial Differential Equations is an important domain of mathematics, so 
semi-classical analysis has also been extremely studied and developed in mathematics. 

1.1 Problematics of quantum chaos 

A common way to express the semi-classical limit or short wave-length limit, is to introduce 
a dimensionless parameter h in the wave equation, called the "Planck constant", which 
corresponds to K ~ l/L where I is the wave length and L the typical size of the domain. 
The semi-classical limit is then H — > 0. 

One possible way to understand semi-classical correspondences is through wave packets 
(or coherent states), which are waves localized as much as possible in phase space, so that 
to mimic the motion of a classical particle. Because of the uncertainty principle AqAp > h, 
in the phase space of position q and momentum p, a wave packet can not be localized better 
than a small domain of surface h, called the Planck cell. This dispersion of the wave packet 
in phase space is responsible for its spreading during time evolution (in a similar way as 
the spreading of a classical disk of surface h). After a finite time (compared to h — > 0), the 
spreading is finite, the wave is still localized. For that reason, one can derive the Van-Vleck 
formula giving the matrix elements of the propagator or the semi-classical trace formula 
[13, 14] [34], and obtain the Weyl formula giving the averaged density of states. 

If the classical dynamics takes place in a compact energy surface, and has exponential 
instabilities (e.g. hyperbolic dynamical systems), the evolution of a wave packet is much 
harder to understand for longer times. This is one of the goals of quantum chaos studies. If 
A denotes the classical Lyapounov exponent of instabilities, an heuristic approach suggests 
that the wave packet which has initial length Aq, Ap > h should spreads into a long and 
thin distribution of length L t > he xt , which rolls up in the energy surface, as soon as L t ^> 
constant t j\og(l/h). This introduces the Ehrenfest time = \ \og{l/h) which 
is an important characteristic time scale in quantum chaos. For longer times, t ^> t^, the 
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distribution may intersect all the other elementary cells of surface h (Planck cells) in many 
distinct components, whose number grows exponentially fast with t. Each component is 
identified with a classical trajectory. As this description suggests (and as been suggested by 
many works in the physical literature [56]), the matrix elements of the propagator operator 
M t between localized wave packets, for t 3> te, could be expressed as a (exponentially 
large) sum over these distinct trajectories. In particular the trace of the propagator could 
be expressed as a sum over Mt oc e xt periodic orbits. This is the content of the Gutzwiller 
trace formula [33] [34]. Although there are numerical calculations that suggest its validity, 
such a semi-classical formula has not yet been proved to be valid for time greater than ^e- 
The difficulty relies in the control of the errors of individual terms which add together. 

A major goal of quantum chaos would be to describe the long time dynamics of wave 
packets up to time tu — l/h (or more), called the Heisenberg time 2 , which would allow 
semi-classical formula to resolve individual energy levels. But at time t ~ tjj the number 
of classical trajectories entering in a semi-classical formula is of order Mt H — exp(X/h), 
extremely large (!), compared to the dimension of the effective Hilbert space which is 1/th. 
One famous observation and conjecture in quantum chaos is that correlations between 
closed energy levels behave as eigenvalues of random matrix ensembles [6] [5]. Many works 
in physical literature assume that semi-classical expressions such as the Gutzwiller formula 
are valid up to time of order tn — 1/^, and deduce some heuristic explanations for the 
random matrix theory, or other important quantum chaos phenomena [22] [35]. For recent 
reviews on the mathematical aspects of quantum chaos, see [17], [18], [59] [62]. 

In this paper we consider a simple model of quantum chaotic dynamics, namely a 
non linear uniform hyperbolic map on the torus, and we show that semi-classical formula 
extends to time scales t = C.tE, where C > is any constant. We do not reach the 
Heisenberg time. The model is particular, but we believe that the methods presented here 
could be extended to more general uniform hyperbolic dynamical systems. In this paper, 
the objective is to discuss the Gutzwiller trace formula, and in a second paper we will 
discuss a second application, the "Van-Vleck formula" which expresses matrix elements of 
the propagator, also for times t = CtE- 

1.2 Characteristic time in quantum evolution of wave packets dis- 
cussed with a numerical example 

In this Section, in order to motivate the importance of the Ehrenfest time, as a qualitative 
frontier in semi-classical evolution problems in quantum chaos, we present and discuss some 
recent results concerning the evolution of quantum states in a hyperbolic flow. One of the 
main challenges in quantum chaos is to deal with both the long time limit t — > oo, and the 
semi-classical limit h — > 0. Usual semi-classical results, such as the Ehrenfest or Egorov 
theorem, concerns h — > first, and t — > oo after. The challenge is to try to reverse this 
order (in order to get information on individual eigenfunctions and eigenvalues), or more 

2 For a time independent dynamics with d degrees of freedom, t H ~ l/h d . For the model developed in 
this paper, tu — 
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modestly, make t depending on h. 

The present discussion is made around the observation of the evolution of a wave packet 
with a numerical example. We refer to Section 2 for a complete definition of the model. 
Let M denotes a non linear hyperbolic map on the torus phase space T 2 , and M the 
corresponding quantized map . A wave packet (a coherent state if) XQ ) is launched at time 
t — 0, at a generic position xo = (<Zo 5 Po) £ T 2 . Figure 1 shows the Husimi distribution (i.e. 
phase space representation) of the evolved state if) (t) = M tr ip Xo at different time t EN. We 
remind that the Husimi distribution of the initial state if) XQ has typical width A — Vh, 
(due to the uncertainty principle AqAp ~ h, and the specific choice Aq = Ap = A ~ \f~R). 
During time evolution, the wave packet center is moving, and its distribution spreads due 
to instabilities of the trajectories. Figure 2 summarizes the main effects we discuss below. 
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Figure 1: Phase space (Husimi) distribution of the evolution of a (generic) coherent state 
at initial position xq = (0.4, 0.2), for different times t = 0, 1, 2, . . . and 2nh = 10 -3 . 
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Figure 2: Characteristic times which appear in the semi-classical limit h — > 0, for the 
evolution of an initial coherent state. t E = j log (l/h)is the Ehrenfest time, (very small "in 
practice") compared to the Heisenberg time tn = 1/^- 



Finite time regime with "no dispersion": We first consider a fixed value of t — Cste, 
and h — »■ (of course t can be arbitrary large in principle). The evolved state if) (t) is 
localized at the classical position x(t) = M t x§. In more precise terms, the semi-classical 
measure of if) (t) is a Dirac measure at x (t) = M t x . The evolved state if) (t) spreads but 
its width is A t ~ e A *A ~ y/H, still of order h 1 / 2 [39] [38] [48]. Because t can be chosen 
arbitrary large a priori, the ergodic nature of the dynamics may have importance if x (t) 
follows a dense trajectory for example. Some well known semi-classical results such as 
the semi-classical Egorov theorem [23], or the Schnirelman quantum ergodicity theorem 
[16] [58] use these finite range of time. 

Linear dispersion regime: Some recent and very general results [13] [4] [36] [9] describe 
the evolved quantum state if)(t), in the linear dispersion regime, which means that non 
linear effects on the dispersion of the coherent state are supposed to be negligible with 
respect to the linear effects. Because the first non linear effects correspond to cubic terms 
in the Hamiltonian, this imposes that A 3 <C h, equivalently e xt f}}' 2 <C h 1 ^ 3 , or t ^t E . In 
our numerical example ^t E = 1.2. 

Localized regime: After that time, the coherent state spreads more and more. But 
its width is still of microscopic size if A t <^ 1, i.e. t \t E . In more precise terms, the 
semi-classical measure of if) (t) is still a Dirac measure at x (£) in that range of time. In our 
example \t E = 3.6. At a time around t ~ \t E , the quantum state has size of order 1, and 
can be described "Lagrangian W.K.B state" [50]. 

Equidistribution regime: For time t larger than \t E , the wave packet spreads and 
wraps around the torus phase space, along unstable manifolds, like a classical probability 
measure. Thanks to classical mixing, a smooth classical probability distribution is known 
to converge towards the uniform Liouville measure for large time. The Husimi distribution 
is expected to behave like a classical measure, and equidistributes, if the different branches 
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do not "interfere" with each other on phase space. After the time ~tE, we evaluate that the 
distance between consecutive branches get smaller and smaller like d ~ e _A (* - * B ' 2 )iintil the 
critical value h is obtained at time t = id,E + tE — I^e (This is indeed the ultimate value, 
because if d ^> h, one can still insert a (squeezed) localized wave packet between two con- 
secutive branches, which means that the branches do not yet interfere). Correspondingly 
to this description, in [7], the authors show that for the linear map, the semi-classical mea- 
sure if){t) converges towards the Liouville measure, in the range of time \tE <C t <C fig. 
J.M. Bouclet and S. De Bievre in [8] obtain a similar result for a non linear hyperbolic 
map, but for t <C \t E . S. Nonnenmacher in [50] reaches the time t <C In [53], R. 

Schubert has described evolution of an initial Lagrangian state under a hyperbolic flow. 
He has obtained similar results, namely equidistribution up to time t <C This is indeed 
similar, because an initial coherent state becomes a Lagrangian state at time \tE- This 
range of time is also considered and controlled in [2, 1]. 

Longer time and interference effects: For longer time very little is known. Some 
arguments and numerical observations in [56] suggest that semi-classical formula applies for 
longer time. This is the subject of the present work and the second paper to come [24]. We 
show that in the range of time t e [0, Cie], where C is any constant, the evolved state if) (t) 
described by its Husimi distribution Hus (x) = \ (x\if) (t)) | 2 , or by its Bargmann distribution 
(x\if) (£)), can be expressed in general as a (finite) sum over different classical trajectories 
starting from the vicinity of the initial state Xq, and ending at time t in the vicinity of 
the point x (similarly to the semi-classical Van-Vleck formula). These trajectories give 
unavoidable interferences effects for time t > \tE-, because one can estimate that the sum 

involves more than one trajectory. Similarly the trace formula expresses Tr ^M*j as a 

sum over periodic trajectories of period t. It is known that in specific cases, revival may 
occur 3 at time t ~ 2t E [27], however it is expected that at least generically, these different 
contributions are somehow uncorrelated, and as a result, the state if) (t) is "generically" 
equidistributed over phase space as can be observed on figure 1. 

An important characteristic time which is not considered here, because far much larger 
than the actual semi-classical approach could reach, is the Heisenberg time t H = 1/h 
(= 1000 in our example). This time is related with the mean separation between eigenvalues 
of M. Some important effect of quantum chaos are numerically observed at this range of 
time, and explained by a Random Matrix Theory approach [5]. It allows to describe 
statistical properties of individual eigenfunctions and eigenvalues. Note that contrary to 
the mathematical works which are "stopped" by the Ehrenfest time, the Heisenberg time is 
extremely discussed in the physical literature, essentially with the random matrix theory. 

3 For a linear map, these interferences effects are responsible for exact revival of quantum states at time 
t = 2tE, and existence of non ergodic invariant semi-classical measure (strong scars), as shown in [27, 26]. 
Reaching the time He for the semi-classical description of a non linear hyperbolic map, was the main 
motivation of this work, although we have not yet any precise idea if strong scarring effect may exist for 
general non linear hyperbolic systems. 



7 



1.3 Results and organisation of the paper 

The main result of this paper is Theorem 12 page 33, which shows that Gutzwiller semi- 
classical trace formula is valid for long time t ~ Clog (l/h), with any C > 0. This formula 



period t, up to an arbitrary small error: for any K, there exists Dc : k > 0, such that 



Tsemi^J = l Det ( D * Mt ~ X ) I V2 eX P i 1 + HE hx,t + ^ E 2,x,t + ■■■ h J Ej^ t ) 



with J > 2 (K + C), and where A^t = j> p<ig — ifcfa is the action of the periodic orbit, and 
Ej >x>t depends on semi-classical Normal forms coefficients of the periodic orbit. 

This formula has already been derived in a similar form (i.e. with Normal form invari- 
ants) many time [32] [61] [43] [55], but for finite time t (with respect to h) giving: 



The main difficulty to extend this last formula for larger time was twice: how to control 
that neither the error term C^k, nor the normal form coefficients Ej )X)t do diverge, when 
t ~ Clog (1/h) is so large. Such divergences could a priori be expected, due to exponential 
increase of the complexity of the dynamics with t. For a linear hyperbolic map, there 
are no semiclassical corrections, the trace formula is exact, and has been derived by J. P. 
Keating [47]. 

In Section 2, we define the dynamics and recall important properties of uniform hyper- 
bolicity. In Section 3, we describe the periodic points and prepare the Trace Formula. In 
Section 4, we show how to control the semi-classical dynamics over large time, which will 
give us a control of the above error term C t ,K- In Section 5, we give a global semi-classical 
Normal form description of the hyperbolic dynamics, derived from a work of David DeLatte 
[19] in the classical case, and as a result, this gives us a control of the above coefficients 
Ej X t uniformly over x and t. Considering the results together, we deduce the semi-classical 
Trace Formula in Section 6. The appendices, that can be skipped in a first reading, give 
proofs of the Propositions used in this paper. The results are illustrated with numerical 
examples in Section 6.2. It is observed there that our bounds on the errors are not sharp 
at all, and that the validity of semi-classical Gutzwiller formula seems to be much better, 
as already suggested in [56]. We discuss again this fact in the conclusion, where we suggest 
an alternative approach with prequantum dynamics, and some perspectives. 

Acknowledgement 1. The author gratefully acknowledges Yves Colin de Verdiere, Patrick 
Gerard, San VuNgoc for numerous discussions, and in particular Stephane Nonnenmacher 
for helpful remarks and suggestions. We acknowledge a support by "Agence Nationale de 
la Recherche" under the grant ANR-05-JCJC-0107-01. 
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2 Quantum non linear Anosov map on the torus 



In this section we describe the dynamical system considered in this paper, a non linear 
area preserving map M on the torus T 2 . This map is constructed as a perturbation of a 
linear hyperbolic map M , and is supposed to be uniformly hyperbolic (which is true for 
small enough perturbations because of structural stability theorem). In order to quantize 
the map M in a natural way, we construct it from a Hamiltonian flow. This section recalls 
some well known results about this construction. 



2.1 Classical Dynamics 
2.1.1 The linear map M 

Consider a quadratic Hamiltonian on phase space x = (q,p) G M 2 with symplectic two 
form uj = dq A dp: 

H (q, p) = -aq 2 + -ftp 2 + iqp, (1) 

with coefficients a, (3, 7 G BL The Hamilton equation of motion for the trajectory x(t) 
are dq(t)/dt = dH /dp = jq + (3p and dp{t)/dt = —dH /dq = —aq — jp. We deduce that 
the flow after time 1 is x(l) = M x(0) with the matrix 

M »='(c B)= e *p(-a - 7 ) eSL(2 ' R) ' P) 

i.e. det(Afo) = AD — BC = 1. We assume that 7 2 > a(3 (equivalently Tr (M ) = 
A + D > 2), so that M is a hyperbolic map with two real eigenvalues e ±A ° where A = 
\J^ 2 — a/3 > is called the Lyapounov exponent. The two associated real eigenvectors 
denoted by uq, sq G M 2 , correspond to an unstable and a stable direction for the dynamics. 
Suppose moreover that A, B,C,D G Z, i.e. M G SL (2, Z). Then for any x eR 2 ,n G Z 2 , 

M (x + n) = M (x) + M (n) = M (x) modi 

so M induces a map on the torus phase space T 2 = M 2 /Z 2 , see figure 3. This map is Anosov 
(uniformly hyperbolic), with strong chaotic properties, such as ergodicity and mixing, see 
[45] p. 154. 



2.1.2 Hamiltonian perturbation 



We introduce a non linear perturbation of the previous map. Consider a C°° function on 
the torus Hi : T 2 — > M (i.e. Hi is a periodic function on IR 2 ). Let Mi be the map on M 2 
given by the Hamiltonian flow generated from Hi after time 1. We compose the two maps 
and define 



M = Mi.M 



(3) 



Which induces a map on T 2 also denoted by M. 
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Remarks: 

1. The final dynamics M* : M. 2 — > M 2 , t G Z, might be seen as generated by a time 
dependant Hamiltonian H (x,t) (periodic in time t G K but discontinuous). 

2. There is a useful relation: 

M(x + n) = M(x) + M (n), VxGM 2 ,raGZ 2 . (4) 



Proof. From M Q (x + n) = M (x) + M (n), and Mi (x + n) = Mi (x) + n, one deduces 
M (x + n) = Mi M (x + n) = M 1 (M (x) + M (n)) = M (x) + M (n), because M (n) G 
Z 2 . □ 



2.1.3 Hyperbolic dynamics and structural stability 

Structural stability theorem: The structural stability theorem states that if the per- 
turbation Hi is small enough with respect to the C 2 norm, then M is an hyperbolic map 
on T 2 conjugated to M by a Holder continuous map H ([3], [45] p. 89): 
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Theorem 2. There exists e > 0, such that if \\Hi\\ c2 < e, then 

M = H.Mo.H" 1 (5) 

with H : T 2 — > T 2 Holder continuous, and M is uniformly hyperbolic in the following 
sense. 

Let DM X : T X T 2 — > T M ^T 2 be the tangential map at x G T 2 . For any x G T 2 , there exists 
a frame of tangent vectors {u x ,s x ) G (T X T 2 ), called respectively unstable and stable 
directions, (chosen such that 

u x As x = l (6) 

i.e. they form a symplectic basis. We can also impose that \\u x \\ = \\s x \\, in order to fix 
the choice) and an expansion rate \ x such that: 

DM X (u x ) = e x *u M {x)i DM X (s x ) = e" Al s M (x), (7) 

and more important: u x ,s x ,\ x are Holder continuous functions of x G T 2 . For any 
x G T 2 , 

< Xmin < A^ < Xmaxi With X m in = mmA^, \ m ax = maxA x . (8) 

1ST 2 isT 2 



In all the paper, we will suppose the perturbation to be weak enough so that the 
previous theorem holds. See Figure 4. 

Remarks 

• Of course, for a null perturbation Hi — 0, M — M , then A^ = A for every x. 

• The conjugation Eq.(5) between M and M is very useful for topological purpose. 
For example, it implies that the periodic trajectories of M (and more generally all 
the symbolic dynamics) are in one-to-one correspondence with those of M . As 
a consequence, M has also strong chaotic properties such as topological mixing. 
However some more refined characteristic quantities of the dynamical systems M and 
M are different. For example the Lyapounov coefficients of corresponding trajectories 
of M and M are different, and this is related to the fact that the transformation 
map H is not C , and does not conserve area. To understand the quantum (or semi- 
classical) dynamics of M, the equivalence Eq.(5) is not strong enough, because in 
semiclassical analysis, we need symplectic conjugations (here, H does not conserve 
area and can not be quantized). 

• In [40], it is shown that u x , s x , \ x are C 2 ~ s functions of x, but we will not use it. 

• For practical purpose, we introduce Q x G SL(2, M), x G M 2 , the symplectic matrix 
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which transforms the canonical basis (e q , e p ) of M 2 to the basis (u x , s x ): 



Qx 



[U x ) q [s x ) q 



Q x depends continuously on i 6 T 2 . We can write: 

DM X = Q MX ( exp ^ ° . 

V exp (-A x . 



(9) 



Q* 1 



(10) 



Explicit expression of u x and s x : Let uo G M 2 be any given vector, and denote 
[uo] G P (M 2 ) its direction in the projective space and [DMJ the action of DM X in 
the projective space. Then the direction of u x is given by 



lim 



which converges for almost every uq. Similarly the stable direction is given by 



lim 



DM M \ (X) 



[so 



where Sq is a given vector. 



2.1.4 Example for numerical illustrations 

In this paper we will illustrate the results by numerical calculations in a specific example. 
The linear non perturbed map is: 



M 



2 1 
1 1 



(11) 



with Lyapounov coefficient Ao, given by e Ao = ^fi = 2, 62 . . .. The perturbation Mi 
is generated by the Hamiltonian 

Hi (q, p) = a cos (2iiq) , a = 0.01 

giving a hyperbolic map on T 2 



M d = Mi M 

Figure 4 shows the torus phase space foliated by stable and unstable manifolds. 



(12) 
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Figure 4: Part of Stable and unstable manifolds of the Anosov map M, Eq.(12), passing 
through the origin (0, 0). At every point x G T 2 , u x , s x are vectors tangent to the foliation. 



2.2 Quantum mechanics 

2.2.1 Quantum mechanics on the plane 

The Hilbert space associated to the plane phase space R 2 is H p i ane == L 2 (R) . The Hamil- 
tonian operator H is obtained by Weyl quantization of H , Eq.(l): 

H = Op We yi (H ) = |<f + ^p 2 + | (qp + pq) , 

Where the operators are defined by (qip) (q) = f qp (q) and (pip) (q) = f — ih (^p^j (q), 

with ip G L 2 (R) in suitable domains, and where the "Planck constant" h > has been 
introduced. The Schrodinger equation defines the evolution of p (t) G H p i ane = L 2 (R): 

dip (t) i * 

and generates a unitary evolution operator between time t — — ► 1, written M : 

y?(l) = M <^(0), M = expU^oJ. (13) 

Similarly, we define ffi = Opvye^z (#i) to be the Weyl quantization 4 of Hi, and Mi = 
exp ^—~Hij the evolution operator after time 1. Finally the quantum map is 

M = Mi. Mq ■ n plane -> H pZ(me (15) 



4 Explicitely, if Hi is expanded in Fourier series: 



Hi (q,p) = c n e a *<™+ n "') 

n=(ni ,n 2 )£Z 2 
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also written M v \ ane in the following, and not to be confused with the quantum map M tor us 
for the torus phase space introduced below. 



Unitary translation operators: For v = (^1,^2) G 

t>2 , irp2 



let 



T 



x 



T v (x) — x + v 



(16) 



be the translation on classical phase space. This translation can be expressed as the flow of 
the Hamiltonian function / (q,p) = (v%p — t^g) after time 1. The corresponding quantum 
translation operator is defined by: 



m def 

l v = exp 



-- (vip - v 2 q) 



(17) 



These quantum translations satisfy the algebraic relation of the Weyl-Heisenberg group[28][51]: 

e iS/h f v+v ,, (18) 



with S — o (v2v[ — viv'2) = —\v A v'. For any matrix M e SL(2, M), one trivially has 
M (x + v) = M (x) + M (v) which rewrites M T V = T Mov M . This intertwining persists 
at the quantum level: 

MoT, = f MoV M . (19) 

We deduce a relation for the non linear quantum map M analogous to the classical 
relation Eq.(4) : 

Mf n = f Mo(n) M, Vn G Z 2 (20) 



Proof. From Eq.(14) and Eq.(16), we can write Hi — ^ 



one has 



-l-nhni T-n' 



0. We deduce that 



0. 



c n T-2nhn- FromEq.(18), 



T n ,,Mi 



0, for any n' G Z 2 



which reflects the periodicity of the function Hi at the quantum level. Then using Eq.(19), 

Mf n = MiM f n = Mif Mo(n) M = f Mo(n) M. □ 



2.2.2 Quantum mechanics on the torus 

At the classical level, the torus phase space was obtained by introducing periodicity on M? 
with respect to the I? lattice, generated by translations T^ o) and 1(0,1). The same con- 
struction can be done in quantum mechanics. The difference is that we have to check that 

(with C- n = cvT, so that Hi is real valued), then 

H 1 =Op Weyl (H 1 )= c„e 42 ^^+« 2 P) (14) 

n—(ni ,ri2)£^ 2 
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the quantum translation operators T(i,o) and T(o,i) commute before we consider their com- 



mon eigenspaces. From Eq.(18), we have T(i,o)T(o,i) 
if and only if h is such that 



"^7(0,1)1(1,0), so T(i, ) , f(o,i) 







N 



def 



2nh 



G N* 



(21) 



We will suppose condition Eq.(21) from now on. We define 



7-Ltorus = jy? G S' (R) / T ( i i0 )V? = <£, ^(o,!)^ = 



In order to have a concrete expression of a wave function <p G Htorus 5 , remark that 
using h- Fourier- Transform defined by (p (p) = f ^7= J dqtp (q) e~ tpq / h , then </? G Htorus is 
characterized by the periodicity conditions (p (p + 1) = <j5 (p) and </? (q + 1) = (q). The 
periodicity of implies that (p (q) = Yln& (? ~~ 77) > with a n G C. The periodicity of 

(g) implies that a n +JV = a n for any n. So G Htorus is specified by {a n ) n=1 ^ N 
deduce that: 



We 



dim7i 



torus 



N 



2nh 



For simplicity we also assume, that N is even, so that, using Eq.(18), T n= („ 1>n2 ) = 
^(1 0)^(01) f° r an y n E 1? . We define a onto projector V : 7i. p i a ne - > H-torus (with a dense 
domain which contains *S (M)), which makes a quantum state periodic with respect to the 
lattice I? in phase space: 



V 



del 



(ni,n 2 )eZ 2 



J (l,0) J (0,l) 



(22) 



From Eq.(19), we deduce: 



MV = VM 

In other words we have a commutative diagram: 

Hplane M 7~~(-plane 

IV IV 

Htorus M T~~Ltorus 

Which means that M induces an endomorphism: 



(23) 



M 



torus 



Mi. M : H 



torus 



H 



tor us 



(24) 



When no confusing is possible, we will write M for this operator M tor us- 



5 7~Ltorus is the well known space for Finite Fourier Transform. 
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2.2.3 Standard coherent states 

We will use in many instances some particular quantum states, the standard coherent 
states [51][57], which are semi-classically localized on phase space. For x = (q,p) G IR 2 , the 

coherent state <p x G H p i ane is f x (</) d = j^tji exp (-z'ff) exp (i*jf) exp f- 4£=|L j • Semi- 
classical localization of <p x at x G IR 2 in phase space comes from the fact that its modulus 
is a Gaussian "localized around" q with width Aq ~ y/h (which vanishes for H — > 0). Its 

ft-Fourier Transform is (p x {p') = ^i/± ex P (~^ft) ex P ex P ^~ ^v^H) 2 } ' s ™^ ar ^ 

localized around p. More algebraically, <p G 'Hpiane (with x = 0), is the ground state of 
the Harmonic Oscillator and is defined by aip = 0, with a — (q + ip) / \/2h. The coherent 
state <p x is then obtained by translation 

ip x = f x ip , x = (q,p) G M 2 . (25) 

For short, we will also write \x) = f ip x for a coherent state. The Husimi distribution of 
a quantum state ip G Ti. p i ane is the positive measure on phase space: 

Hus^ (x) d = — |(a#>| 2 



The closure relation is ([51] p. 1 



Id /W , = / — t\x){x\ (26) 
A coherent state on the torus is defined by periodicity, using Eq. (22): 

| •£/ torus 'P\x S j G T~Ltorus 

They provide the closure relation: 

f dx 

Id/w torus = / ^-r\x) toms (x\ (27) 
i[o,i] 2 2vrn 



3 The map as a sum over periodic points 

In the first part of this section we characterize and count the number of periodic points of 
the map M with a given period t G Z. These points are shown to play an important role 
in the second part, where we will express M\ orus , defined in Eq.(24), and its trace in terms 
of them. Some parts of this section can be found in [46] or [52]. 
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3.1 Periodic points of M 

Consider a discrete time t G Z\ {0}. A point x G IR 2 gives a periodic point [x] G T 2 of M 
with period t if 

M*[x] = [x] 

^> 3n G Z 2 / M*x = x + n 
& 3n G Z 2 /T t (x) = n 



with the map 



T A — 

it — 



R 2 -> IR 2 
x — >• T t (x) = M 1 (x) - x 



Periodicity of %\ From Eq.(4), we get 

% (x + m) = % (x) + T m (m) , Vx G IR 2 , Vm G Z 2 



Where 7o,(t) (x) = f Mq (x) — x. In particular 7o,(t) (m) G Z 2 for any m G Z 2 . We introduce 
the sub-lattice of Z 2 : 

A, T m (Z 2 ) C Z 2 (28) 

generated by e\ = %^ t ) (1,0) G Z 2 and e 2 = 7o,(t) (0, 1) G Z 2 . We denote by C« C Z 2 the 
elements which belong to the origin cell of A t : 

C/= f {nGZ 2 / ^(n)e[0,lf} 

Of course C t = Z 2 /A t . 



Proposition 3. For a small enough perturbation, % is a diffeomorphism, for every t. 
Periodic points with period t can be labeled by n G Z 2 : 

x n , t = Tf 1 (n), ne Z 2 (29) 

Different values n, n' may give the same point \%n,t\ — [-^n'.t]- The periodic points [x n ^] G 
T 2 on the torus are in one to one correspondence with n G C t . The number of periodic 
points with period t is 

M = tf (C t ) = |det (T , (t) ) | = |det (Ml - I) | (30) 
e Ao '-2 + e" Ao ' (~ e v /ort»l) 
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Proof. Using Eq.(lO), the matrix of (DT t ) x is 

(D%), = (DM% - U = Q M „ ( < A -» exp ( ° ^ ) - H 
with A^t = £)£l -^Af*'x > so ^mint < A Z) t < A max t. We have supposed that X min > 0. Then 
det ({DT t ) J = det ( ( e *P <^> ^ ( ° ^ ) - fl^fl.) 

The matrix Q~}Q X goes to identity (uniformly in x, x' G T 2 ) when the perturbation van- 
ishes. Therefore, we can write Q~}Q X = Id + eA x>x r, where A XtX > has matrix elements 
bounded by 1 in absolute value, and e goes to zero when the perturbation vanishes. One 
computes 

det {{D%) x ) = 2 (1 - cosh A*, t ) + eA n (e~ A ^ - l) + sA 22 (e A *<* - l) + e 2 det (A) 
and deduces 

| det ((DT t ) x )\ > 2 (1 - e) (cosh(A min t) - 1) - 3e > 
for small enough e. We deduce that % is a diffeomorphism onl 2 . □ 

Remarks: 

• Note that two periodic points x n ' jt , x n<t may belong to the same periodic orbit of 
M torus . Explicitely: 

3V, t = M (x n>t ) <&ri = M (n) 

so periodic orbits of M on T 2 with period t, are in one to one correspondence with 
periodic orbits of M acting on the finite set C t = I? / K t . 

• Here is the number of periodic points with period t, for the example Eq. (11), see 



figure 5: 
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3.2 Expression of the quantum map M^ orus , using periodic points 

For a given t G Z, we consider 



plane J • ri plane iltorus 

defined on a suitable dense domain (which contains the Schwartz space S (K)), where V is 
defined in Eq.(22), and M p i a ne is defined in Eq.(15) . This operator [PMpi ane ) is important 
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Figure 5: Periodic points of map M, Eq. (12), with period t = 7. There are Aft = 841 
periodic points. 



to look at, because if one wants any matrix element of the quantum map (il)2\Mj: orus \i/ji), 
with |^x), \i/) 2 ) £ 7~Ltorus-, then one just has to consider "lifted states on the plane" 
1^2) G Hpiane, such that = V\ipi), i = l,2, and then 



We will now write M instead of M p i ane . One writes: 




with 



: U 



plane 



n 



plane 



where n G Z 2 , and t G Z. The corresponding classical map is of course 

i£ 2 -> M 2 
x -> M n f (x) = M* (x) - n 



M 



(31) 



(32) 



(33) 



(34) 



The map M n >t will be used many times in this paper. It is an hyperbolic map (translation 
of M t ), whose unique fixed point is the periodic point x„ )t , given by Eq.(29), because: 
M nyt (x n>t ) = M t (x n ,t) -n = x n>t . 



3.2.1 Periodicity of the decomposition 

We decompose now the sum over n G I? in Eq. (32), with respect to the lattice A t defined 
by Eq. (28). First, any n' G Z 2 can be decomposed in the unique way: 

n = n + 7q (m) , n G Ct, m G 7? 
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Observe that for n, m e I? (we use Eq.(20)), 

M n+ T (m),t = T- n -T a {m)M = T_ n _jv/ ( m ) +m M* 

= e^^ m )/ n f m f_ n f_ Mo(m) M* = e^ m V R f m f_ n M t f_ 

%T(n,m)lhrf jTr rp 
° ± m lvl n,t- L — m 

The phase T comes from Eq.(18) and is given by: 



T (n, m) = — (n A Mq (m) — m A (n + Mo (wi))) 

But 2JF is an integer so e i:F ( n < m )/ h = e ii-xNT{n,m) _ ^ because we have supposed N even. 
Therefore 

(35) 



with 




M Ct M n , t 

n&Ct 



From Proposition 3, this last expression of Mc t is a finite sum over periodic points, with 
Mt terms. 

3.2.2 Trace of 

Remark that Tr ^M* orus J is well defined because M t ' ortts acts in Titorus which is a finite 

dimensional space. On the opposite M n t = T^ n M t is a unitary operator in TL p i an e = L 2 (M) 
and is not of Trace class. Nevertheless we will define a linear functional which can be 
thought as a trace. 



Proposition 4. For any t <E Z, and any n G Z 2 ; the following integral is absolutely 
convergent: 

def 



TIM, 



n.t 



dx 
2Tih 



(x\M nit \x) 



(36) 



where \x) is a standard coherent state at x el 2 defined in Eg. (25). We have the relation 



Tr M 



torus 



(37) 



neCt 



Eq.(37) is an exact formula (not semi-classical). It is a sum over M t terms. This formula 
does not use the assumption that M is hyperbolic. 
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Proof. A coherent state \x) belongs to S (R). Using closure relations Eq.(26) and Eq.(27), 
together with Eq.(31) and Eq.(35), we compute: 

Tr M* oru J = / ( x \Ml rus \x) torus = / -—(x|PM* Zane |x) 

^ ^ nGCt mGZ 2 

= Y] [ y^ ^-(x-m\M n!t \x-m) 
\n 11 2 27Tft 

The sums are absolutely convergent. In particular: 



1°'^ mez 2 

and we obtain Tr ( M* orus ) = £ ngC , T ( M n „ t ) . 



4 Semi-classical description of the dynamics 

We do not have yet considered the semi-classical limit h — > 0. In this section, we give the 
essential Theorem (Theorem 6 below) which shows that in the semi-classical limit, long 
time quantum dynamics can be expressed in terms of individual classical trajectories. 

4.1 Neighborhood of a point in phase space and localized states 

We first introduce the notion of a semi-classical neighborhood of a point in phase space. 
Let x G R 2 be a point in phase space, and < a < 1/2. Let 

D x , a = {yeR 2 /\y-x\<h l ' 2 - a } (38) 

be the disk of center x and radius fr 1//2_Q , which shrinks to zero as h — > 0. We define the 
truncation operator: 



Px,a = [ Pr\y){v\ = OpAW (XD x , a ) 

JytD x , a 



(39) 



being the Anti-Wick quantization of the characteristic function of the disk D xa [51]. We 
will often drop the index a, and write P x = f P Xja - In the case x = 0, we will drop the 
index x, and write: P a == P x =o,a- Notice that in comparison with Eq.(26), the integral is 
restricted to the disk D x a . The meaning of the operator P x is that it truncates a quantum 
states in the vicinity of the point x. A quantum state is said to be localized at point x, if 
this truncation has no effect on it in the semi-classical limit. More precisely: 
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Definition 5. 

depends on h - 



A sequence of normalized quantum state iph G Tlpiane (sequence which 
> 0) is said to be a— localized at point x G M 2 , if 



L 2 



o(h°°) 



(40) 



Examples: a coherent state \x) is a— localized at x, for any < a < 1/2. For any fixed 
n G N, let \n) be the eigenstate of the harmonic oscillator: \ {p 2 + q 2 ) \n) — H (n + |) \n), 
and let \x,n) = T x \n). Then |x, n) is a— localized at x, for any < a < 1/2. 



4.2 Semi-classical evolution in a neighborhood of a classical tra- 
jectory 

Let x G M 2 , and consider the classical trajectory x (t) = M f x G M 2 , for t G N. The following 
Theorem will be useful in order to compute matrix elements of the quantum evolution 
operator, of the form (^ 2 |M* \ipi), where ip\ is localized at x (0), and ip2 is localized at x {t). 
From Eq.(40), (^|M*|^i> = (^| A,(t)M*A-(o) l^i) + 0(h°°), so the computation involves 
the operator P x{ t)M t P x{ Q ) . 



Theorem 6. For any C > 0, any t, such that 1 < t < Clog (l/h), any < a < 1/2. any 
K > 0, there exists Ck > 0, such that for any x = x (0) G M 2 , then in L? operator norm: 



P x ®M t Px(o) ~ P x{t )MP x {t-i)MPx(t~2) 



MP, 



x(0) 



L- 



< c K h K 



(41) 



The proof is given in appendix A page 37, but we give the main idea below. 



Remarks: 

• To shorten we write that the error is O (h°°) uniformly with respect to x G M 2 . 

• This relation means that in order to compute P x ^ t )M t P x ^, we can introduce trun- 
cation operators all along the trajectory, without changing the result significantly in 
the semi-classical limit. In other words the operator P x ^ t )M t P x ^ depends only on the 
dynamics in the vicinity of the trajectory x (0) — * x (t). This result will allow us to 
use normal forms in the next Section, which give a nice description of the dynamics 
in the vicinity of any trajectory. 

• Idea of the proof: Let us explain here the main idea of the proof but at the level 
of classical dynamics. The proof follows the same idea. From definition Eq.(39), the 
operator P x ^ = f xeD t^\x){x\ truncates quantum states on the disk D t of small 
radius /I 1 / 2- ", and center x (t). The transcription of Eq.(41) in classical dynamics is 

A n M* (D ) = D t n M (A-i • • • M (A n M (A))) (42) 
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To show this, let G\ = M (A) \ A- Let F 1 such that M 2 = A U F 1 U G x is a disjoint 
union. So A n M (A) = 0. Then 

A n M* (A) = A n M* _1 (M 2 n M (d )) (43) 
= An M t_1 ((A u Fi u d) n M (A)) 
= (A n M*- 1 (A n M (A))) u (A n M* -1 (GO) 

Now observe that Gi is the set of points coming from A but who leave the set A in 
the unstable direction. Due to uniform hyperbolicity and the fact that we consider 
the dynamics on the cover M 2 , the set G\ goes away from the trajectory x (t) in 
the unstable direction. This gives: A H M l ~ l (GO = 0. Then Eq.(43) simplifies 
to A H M* (A) = A H AA 1 (A n M (A))- Repeating this argument, we get Eq. 
(42). (For illustration, see Figure 9 page 41). 

• From the idea of the proof given above, it is clear that Eq.(41) holds because it 
concerns the phase space cover IR 2 . Points who leave the trajectory in the unstable 
direction never come back. At the semi-classical level, there is no interference effects. 
The interference effects come when passing to the torus which is compact, and are 
due to the sum Eq.(32). 

Consequence for the trace There is a consequence of Theorem 6, which shows that 
the integral Eq.(36) up to a negligible error, can be restricted to a neighborhood of the 
fixed point x n)t of the map M n j , Eq.(34). This Lemma will be useful in order to obtain 
the semiclassical Trace formula. 

Lemma 1. For any C > 0, any t, such that 1 < t < Clog (l/h), any < a < 1/2: 

T (M M ) = / pz(x\M n>t \x) + O (h°°) (44) 

= Tr (MnA n ,) + O (fi 00 ) = Tr (P Xn , t M n ,P Xn .) + O (h°°) 

The error is uniform with respect to the point x n ^ (i.e. with respect to t and n G Z 2 ). 

Note that P x is a trace class operator, M n4 is bounded, so (M n;t P Xn ^j is also trace 
class. Although intuitive, our proof of Lemma 1 is quite long and is given in Appendix B. 



5 Semi-classical non-stationary Normal Form 

In the previous Section, we have explained that long time quantum dynamics can be ex- 
pressed with the operator A( t )M'A( ), where x (t) = M l x (0) is a given classical trajectory. 
Using truncation operators all along the trajectory, we have shown in Theorem 6, that the 
operator A^M* P x {o) depends in fact only of the vicinity of the trajectory x (t'), t' G [0, t}. 
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There appeared operators like Pm( x )MP x , with x G M. 2 . This suggests to use a local de- 
scription of the dynamics along the trajectory, in terms of a Taylor expansion up to a 
given order J. Using convenient canonical coordinates, we can make this description in its 
simplest form, called a normal form expansion. This is Theorem 8 below. A normal form 
expression of the dynamics is particularly interesting for large time, because along a given 
trajectory we will obtain a product of normal forms operators which is quite easy to calcu- 
late (because they commute together). In particular, for a periodic orbit the computation 
of the trace will be explicit. 

For classical hyperbolic map on the torus, David DeLatte in [19], has shown that there 
exists a global and uniform normal form expression, called non-stationary normal form, 
which is unique up to co-boundary terms. In this section we develop the semi-classical 
version of his result, and use it to control long time evolution. Semi-classical normal 
form for individual unstable trajectories is already a well known and very useful tool in 
semiclassical analysis, see [11, 15, 12], [60, 61], [54], [42, 41]. Our approach of semi-classical 
non stationary normal forms follows closely the exposition of J. Sjostrand in [54], but will 
be adapted to the normal form approach of David DeLatte [19], which is uniform over 
phase space and not individual for periodic orbits. Therefore it will give a satisfactory 
control of the normal form expressions for long (periodic) orbits. 

In this Section we give the main result useful for our purpose, and in appendix D, 
we give the proofs and more details, in particular we present there the semi-classical non 
stationary normal form theory in terms of Hamiltonian flows. This could be useful in order 
to extend the present results to more general hyperbolic Hamiltonian flows. 

5.1 Semi-classical non-stationary normal form 

We give here the semi-classical version of a Theorem of David DeLatte [19] about non 
stationary normal forms. 
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Theorem 7. Let J > 2 even, be the order of the normal form calculation. Let < 
a < 1/2. The evolution operator M in the neighborhood of any point x G M 2 ts well 
approximated by a normal form operator N x : 



Pm( x) mP x - P M{X) [t M{x) N x f- x \ A 



< Ch A (451 



with A = — a') ( J + 1) — 1, any a < a' < 1/2, and C independent of x. The operator 



N x = exp {-J**) ( 46 ) 

is generated by a Hamiltonian with a total Weyl symbol K x (g, p) which is a normal form 
up to order J: 

K x (q,p)= Yl X U,(*) hl (WY ( 47 ) 

0<l+j<J/2 

-Vi.O) = is l° ca l expansion rate already introduced in Eq.(7), and the meaning of 
the other Xij,( x ) is discussed below. 7~ x is a product of unitary operators: 

% = f x Q x U G ^U G2x . . . U Gnmax ^ (48) 

where T x is the translation operator Eq. ( 1 7), Q x is the Weyl quantization of the linear 
symplectic map Q x , Eq.(9). The next operators give non linear corrections: the operator 

?7g„,x = ex P (~ jpn^ is generated by G n , x whose Weyl symbol is equal to 

G n ,x (q,p) = gi, a ,b,(x)h l q a p h , with 3<2l + a + b<J 

in a neighborhood of the origin (the index n enumerates all the indices (/, a, b) in the range 
3<2l + a + b<J, and with increasing order of {21 + a + b) ). The functions gi, a ,b,( x ) an d 
\,3,{x) f or (l>j) 7^ (0,0), are continuous with respect to x G T 2 . The function Ao,o,(a;) ^ s 
continuous with respect to x G M 2 (but not periodic). 



The proof is given in appendix D. 



Remarks: 

• Eq.(45) is interesting if the error vanishes in the semi-classical limit K — > 0, so if 
(J + 1) (| - a) - 1 > & J > for example if J = 2 and < a < 1/6. 

• If J = 2, the normal form description is just at the level of linear approximation. It 
is the quantum version of Eq.(10) and relies on uniform hyperbolicity. In the proof, 
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it will be clear that the next order terms are obtained iteratively from this linear 
approximation, as usual in normal form calculations. 



• Eq.(45) has a classical version, with maps N x , T x instead of operators N x , T x re- 
spectively, and with similar Taylor expansions, but without the semiclassical terms 
h l , I > 1. The classical normal form is expressed in Figure 6, and corresponds to the 
theorem 1.1, p. 238 given in [19] by David DeLatte (the formal version). In [19][20], 
David DeLatte proved a more stronger result: if M is analytic, then normal form 
series of K x and % are convergent in a neighborhood of (q,p) = (0,0), for J — > oo. 
Thanks to truncation operators, we will not need this result. 




Figure 6: This picture traduces the conjugation relation Eq.(45) of the non stationary 
normal form. Here, in a neighborhood of a point x, the classical map M = T M ^N x T~ l 
is conjugated to a map N x which is a normal form up to order J, and has as hyperbolic 
fixed point. 

• The main non trivial result in Theorem 7 is the continuity and periodicity of N x and 
T x with respect to x. This is a global constraint over the torus, and relies on the 
uniform hyperbolicity of the dynamics. This continuity is already given in Theorem 
2 for the expansion rate A^. In the proof, this appears in Lemma 12. 

• -Vo,(aO 1S ca U e d the Action of the path x — > M(x). It is given by the integral 
•Vo,0) — — / \ (pdq — qdp) — Hdt along the trajectory, as explained in Eq.(105) page 
60. Ao,o,(a;) is a constant term in the function Eq.(47) and does not appear in the 
classical version of the Theorem, but is fundamental to explain the interference effects 
in quantum mechanics. For a geometric interpretation of Ao,o (x) in terms of parallel 
transport on a Complex line bundle over the phase space (called the prequantum 
bundle), see [25] and references therein. 

• We will use later on, the fact that for (J, I) ^ (0,0), \j,i,( x ) is a continuous function 
of x G T 2 , and is therefore bounded: 

Aj,Z,mm = miOrgT 2 (^j,l,(x)) < ^j,l,(x) < \,l,max = max xgT 2 (A^^)) (49) 
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The function A^o,^) = A x . is equal to the expansion rate introduced in Eq.(7), and 
is called the Lyapounov cocycle. The function A 2j o,(x) of x, is called the Anosov 
cocycle and is an obstruction to C 2 regularity of the unstable/stable foliation, see 
[40], or [37] p. 289. These two cocycles are the first terms of the series \jj of 
(semi-classical) cocycles. 



5.2 Semi-classical normal form for a long orbit 

We will now combine Eq.(41) and Eq.(45) along a given trajectory x (t) = M l x starting 
from i6l 2 . Let us first define 

def tC 

N Xjt = N M t-i( x )N M t~2( x ) . . . N x 

to be the product of Normal forms N x = exp ^— \K X J along the trajectory. Quantum 

normal forms Hamiltonian K x at different point x commute together (this is proved in 
Eq.(67) page 31). So the product N Xit can be written 



N X)t = exp ( --tK X)t 



(50) 



where K Xjt = Opweyi (Kx,t) has total Weyl symbol: 

dcf 1 

Kx,t(q,p) = - ^Kmsx^p) 



which is just a time average along the trajectory. Using Eq.(47) we can write 



Kx,t(q,p)= Yl Vij,(x),th l (qpY 

0<Z+j<J/2 



(51) 



with coefficients 




(52) 



From Eq.(49), we deduce that [AjjMj is bounded uniformly with respect to i 6 T 2 and 
t G Z (i.e. for any orbit): 



A 



j,l,min 



(53) 



We can now state: 
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Theorem 8. For any J > 2, any C > 0, any < a < 1/2, there exists C\ > 0, such that 
for any < t < C\og(l/h), and any initial point x G M 2 , any 1/2 > a' > a, 



PM t {x)M t P x — Pm*(x) V^M^x^x^x 1 ) Px 



< t 



IP- 



(54) 



As a corollary, for a periodic orbit, i.e. any fixed point x n j of Ml orus given by Eq.(29): 



(55) 



where M n j is defined in Eq.(33), T[M n A is defined in Eq.(36), and 



a def 1 
Ai,t = T. n A X n, 



(56) 



Remarks: 



The constant A n j and jUo,o, ten t),t = ~ j f \ (pdq — qdp) — Hdt defined in Eq.(52), 
contribute to the total action of the periodic orbit defined by: 



-4-n,t — A n)t — t/io,0,(a;„,t),t 



(57) 



Proof. The proof of Eq. (54) combines Theorem 6 and Theorem 7 in three steps. First from 
Eq.(45), we deduce that 



P x ®MP x{t -x)MP x{t - 2) . . . MP m 



Px(t) [%(t)N x (t-l)'T xi l_ 1) ) P X (t-l) [%(t-l)N x{t -2)% { l_ 2 )) Px(t-2) ■ --PxiO) 



(t-2) 



I? 



GtW 



We can now apply Theorem 6 to the sequence of hyperbolic maps (Tm( x )N x T x x ) in order 
to take off the intermediate operators P x ( s )- We obtain: 

P x (t) {^M t (x)^x,t^ x P x (o) 

-Px(t) (%(t)N x{ t-l)T- i t_ 1) ) P x (t-1) \J X {t-l)N x (t- 2 )f~ { t_ 2) ) P x (t-2) ■ ■ ■ Px(0) L2=0 (h°°) 



Combining the last two equations with Eq.(41), we obtain finally Eq.(54). Now we will 
prove Eq.(55). First Eq.(54) for x = x n>t gives 



P M\x n , t ) M ' ' P xn,t ~ PMt(x n , t ) (yMt(x n ,t) N x n ,t,t?x n \) Px n ,t 



ft 



28 



But M l x n)t = x n j + n, so Af*(x„, t ) = f n P Xnt f_ n . From Eq.(48), and Qm*^ = Qx nt (due to 



periodicity), we have T M t (Xn t) = T MlyXnt) T^ Xnt T Xnt = e lAn ' t/h T n T Xn t , with A n , t = \n Ai„,i. 
The last equality comes from Eq.(18). We obtain: 

\P Xn , t M n ,P Xn , t - e^' h P Xn , t [% nt X n ,A-]) K\ L2 < tC i ^ J+1 )(i- Q ')- 1 
Lemma 10 page 48 allows us to pass from operator norm to Trace norm estimates. It gives: 

Tr (P Xn ,MnA n ,) ~ ^- tA Tr (P Xn , [% n>t N Xn ,A^) K) \ < *d ft< J + 1 X*-«')-i-«' 

(58) 

We want now to take off the operator T x t . Remind that T x is a smooth canonical 
transformation which send point to i. If < (3 < a, we have in operator norm 

P X ,f3%P ,a = P x ,pf x + 0{h°°), P , a %P Xt p = t x P 0> p + O (ft°°) , and P x ,f3P , a = P x ,p + O (h°°). 

So, with (3 > 7 > a, 

Tr (Px n , t ,a {^x n , t N XnU t'f x 2t) Pxn,t,a) = Tr (P*n,t,a (% nit Po,jN Xnut P ^T Xn ^ P XnM <^j + O (ft™ 

(59) 

= Tr [P Xn ^ (f Xnd P Qa N Xn ^n^ Kup) + O (ft* 
= Tr {f Xn ^K,aKA~^ + O (ft 00 ) 
= Tr(p Oi7 iV Xnitit P o , 7 N )+0(^ 00 ) 



For the second line, we have used the property that the trace does not depend on the choice 
of 0, up to O (ft 00 ). Finally, Eq.(58,59,44) together give Eq.(55). □ 

5.3 Post-Normalization 

5.3.1 Post Normalization corrections 

The Weyl symbol K x t (q,p) given in Eq.(51) is a function of the product (qp) only. Let 
us write it K x j (qp) ■ The quantized operator we have to consider in Eq. (50) is K x ^ = 
Opweyi (Kx,t (qp))- First remark that, in order to compute the trace of the propagator, or 
its matrix elements, it is easier to deal with a function of the single operator Opweyi (qp), 
and second, Op Weyi (K X)t (qp)) ^ K x>t (Op Weyi (qp)) in general (see e.g. Eq.(65) below). So 
we have to find a function K xt of a single variable, such that 



K Xit = Op Wey i (K Xtt (qp)) = K X)t (Op We yi (qp)) 



(60) 
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Proposition 9. If K x ^(qp) = Ylo<i+j<j/2l Jj hj,x,t^ 1 {qpY *s expressed as a power series in 
(qp) and h, (as we get in Eq.(51)), then K x>t is given by: 


K x , t {qp) = ^ hhx,th l (qp) 1 

l,j>0 


(61) 


where the coefficients jj>i,j, x ,t are given explicitely from p>i,j, x ,t by: 




fii,j,x,t = E f,mP>l',j>,x,t, I' = I - 2m, f = j + 2m 

0<m<[//2] 


(62) 


where Eji^ m is a numerical factor determined by recursive formula: 




3 2 

Ej+l, m = E j>m — — i£y-l),(m-l)> V j, m > 1 




Ejfl = 1, E 0)Tn = 0, Ei >m = 0, for m > 1 


(63) 



Remarks: 

• In particular, at the principal and sub-principal level (i.e. order I — 0, 1 in h l ), there 
is no change: i.e. jj>ij, x ,t — Hij,x,ti if / = or Z = 1. 

• In [31], paragraph 6, Alfonso Garcia-Saz gives such expressions as a special case of a 
more general problem. See also appendix I in [10]. 

Proof. We use the Moyal start product Jj defined in Eq.(88), which is equivalent to the 
product of operators. In other words, we want to express (qp) 1 in terms of product of 

monomial terms (qpf 3 = f (qp) $ (qp) j} . . . jj (qp) . We just apply Moyal product formula 
Eq.(89), with operators A = qp, B = (qp) 3 , and observe that (AV n B) — 0, for n — 1 and 
n > 3, and (AV 2 B) = -2j 2 (qp) 3 ' 1 . This gives 

(qp) ft (qp) 3 = (qp) 3+1 + h 2 ^ (qp) 3 ~ l (64) 

Or equivalently 

Opweyl (qpf +1 = (Opweyl (qp)) (Op Weyi (qpf^j - (Op Weyl (qp) 3 ^ (65) 

We deduce that 

(qp) 3 = E i,mfr 2m (qpf 3 ~ 2m) (66) 

m=0 
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This last formula says that each resonant term of Opweyi (Kx,t) can be written 

b'/2] 



l^l,j,(x),th l Op Wey l ([qpYj = Hl,j,(x),t ^2 E j,™ft l+2m (Opweyi (qp)Y~ 

m=0 



2m 



as desired. 

Proposition 10. (qp) n and (qp) m commute with respect to the star product: 

[(qp) n , (qp) m \ = (qp) n (t (qp) m - (qp) m H (qpf = 
therefore if F,G £ C°° (R), 

[Opweyi (F (qp)) , Op Wey i (G (qp))] = 



□ 



(67) 



Proof. By induction on the power: suppose that \(qp) n , (qp) m ]» = 0, for n,m < N. Using 
Eq.(64), we compute 



(qp) N+1 , {qpf 



(qp) N+1 tt (qp) 11 - (qp) n tt (sp) 



(sp) tt (sp)* (<?p) n -ti 21 — (qp) 1 " 1 tt (gp) n 



TV 2 



(?p) n tt (ep) tt (sp)* + ^ 21 r- (?p) n tt (?p) 



iV 2 



vJV-l 







from associativity of the star product. There is a more direct proof, using the post- 
normalization transformation F — >■ F described above: 



[Opweyi (F (qp)) , Op We yl (G (qp))} = F (Op We yl (qp)) , G (Op W eyl (qp)) 



0. 



□ 



5.3.2 Semi-classical Trace formula for a periodic orbit 

Semi-classical expansions for the trace of hyperbolic normal forms is explicitly done in 

[43]. It can be used to give Tr { i P a N Xn ,tPa S j in terms of the semi-classical post-normalized 

cocycles jj>ij, x „,t defined in Eq.(62). We do this in appendix E, Proposition 24, page 62. 
We can deduce the following proposition: 

Proposition 11. Let J > 2. For any C > 0, any < a < 1/2, there exists G\ > such 
that, for any time < t < C\og(l/h), any n £ C t , one has a semi-classical expression: 



Tr ( PaN Xn t tP a j T sem i j (x n j 



(68) 
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where x n j is the fixed point of the map M n> t, Eg. (34), and 



T S emi,j (xn,t) = exp ^-a ^ ' ^ ) exp (-zt/ii )0)!BBl t, 



h 

with a semi- classical expansion 

£ Xnt = l + hE l + h 2 E 2 



2 sinh ( 



H0,1, x„,tt\ x ^,t 



([J/2]-l) 



where E s depends on t and jlij jXn! t , with (I + j) < s + 1. We have the bound 

\Es\ ^ t E max s 

where E maxs does not depend on t nor on n. With Eg. (55) we deduce that 



M n ,t) -e^^ emi)J (z n , t )| < dh [J/2] 



(69) 



Remark: Explicit expressions of E S: s > 1 are quite complicated for large s, and are 
given in Eq.(115) page 63, and Eq.(123) page 65. It appears in appendix E that with Weyl 
quantization then H\fl tXn ,t = (and /^j = for / odd). We will use this simplification 
below. 



Proof. Proposition 24, page 62, gives 



Tr ( P a N Xn t iP a j T sem i j [x n t 



2 sinh 



< dB, where d 
E im hwA given 



does not depend on t and n, and with B = max n f ( 

by the maximum over periodic orbits of the next order term in the series of T semii j (x n j). 
The uniform control of the semi-classical cocycles, over the periodic orbits, Eq.(27) gives 

,l,min ^mim 

where X min > is defined in Eq.(8). Thus (2 sinh ' 

-x. 



< 



(j 2e -x min t/2 < Q 2 h\ 2t E ) , because we can write e"^™*/ 2 = e ~ A ° iB V J = fr\ 2*i ; The 
uniform control of the semi-classical cocycles together with the bound Eq.(llO) also gives 
that -%//2] < t s E maXyS where E maXjS does not depend on n, t. Now if \t\ < O (\og(l/h)) 

then \t s \ < tT £ for any e > 0, so at final dB < Q^^) 'h^H^ < dfr [J/2] , if we take 



i I x 



This gives Eq.(68). 



TIM, 



e^ h T semh j(x n , t ) <dt* J/2] 



From Eq.(68) and Eq.(55), we deduce that 
c ^ h {J+i){\-<*')-i-*' _ Wg observe that for large Jf this 

gives a bound 0(h°°). So we de- 
duce that the actual bound is given by the next order term in the series of T semi< j (x n> t), 
namely B = O (fii J / 2 ]) already computed above: in more precise terms, this is because 

T (M n , t ) - e lA -/ h T semiJoo (x n , t )\ = O (h°°), and \T semiyJoo {x n , t ) - T semi ,j {x n , t )\ = O (ti J M). 



□ 
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6 Trace of the quantum map 



6.1 Semi-classical trace formula 

The following Theorem is one of the main result of this paper. It gives a semi-classical 
expression for the trace Tr ^M| orus J , as a finite sum over the fixed points x n>t of the classical 
map M l torus . 



Theorem 12. For any K > , any C > 0, any J > 2 (K + C), there exists C\ > 0, such 
that for any time \t\ < CtE, with tE = j- log |, any admissible value of h (see Eg. (21)), 

the trace of M\ orus is given by the semi- classical formula: 



Tr M' — T. 



torus I 1 semi,t,J 



< C\h K (70) 



w- = E «p (^) j-^ks) < n > 

where each term of the finite sum is associated with a fixed point x n .t of M\ orus given in 
Eq.(29). A nt t = A n —fio t0tXntt t is the action of the periodic orbit, defined in Eq.(57). fJ>i,j, Xn ,t 
and fii,j,x n ,t (ire the semi- classical normal form coefficients (cocycles) of the periodic orbit 
computed up to a given order J > 2 (i.e. with 2(J + l) < J). S Xnt is a semi- classical 
series: 

£ *n,i = 1 + HSE ° = 1 + HEl + H2E * + • • • + ^ [J/2l ~ 1 ^[J/2]-l 

l<s<[J/2]-l 

where E s depends on t and fiij tXn ,t (with (I + j) < s + 1), and in particular is bounded by 

\E I < t s E 

\ J -'s\ — J -'rnax,s 

where E maxs does not depend on t nor on x n)t . 



Proof. We use Eq.(37), which is a sum with J\f t = e 



-A t 



terms (from Eq.(30)). 



We write Af t < e A °* = e ~ XotE{ - t/tE) = h { ~ t/tE) < fr° '. Each term T (M n ,tj is approximated 

by a semiclassical expression given in Eq. (69) , with an uniform error dh^W. By a simple 
triangular inequality (i.e. we sum the magnitudes of the error bounds) we deduce that 



Tr M 



torus 



T f 



emi,t,J 



This is O (h K ) if J > 2 (K + C) (for J even). 



□ 
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Remarks: 



• If we look at the whole proof, the limitation to these logarithmic time (any multiple 
of the Ehrenfest time t E = ^ log |) is present many times. But the main limitation is 
due to exponential proliferation of periodic orbits in the hyperbolic system, At — e A °*. 

• The bound on the errors could be improved at many places in the proof, so the 
condition J > 2 (K + C) is not sharp. But the main limitation in our proof is that 
we simply estimate the total error as the sum of the absolute values of the error from 
each periodic orbit. However its seems that all these errors compensate each other, 
as observed in the next Section. However, a rigorous analysis of these compensations 
is out of reach for the moment. 

• Let us comment again on the estimation used in the proof. We have (2 sinh ( w,1 '^ n,t )) 
g-Aot/ 2 f or i ar g e ^ S0) ^ we wr jte Eq.(71) as T semi ^ y j = J2 n ec t ^ni, each term associ- 
ated to an individual periodic orbit |T n t | < e~ A() */ 2 decreases, but Mt — % (C t ) — e A °* 
increases faster. We get the bound \T aemi ^j\ < e Xot ^ 2 which is greater than dim (7-^) = 
N = 1/ (2ttH) for t > ^ log (l/h) = 2t E . This shows that for t > 2t E , there are nec- 
essarily some cancellations among the complex amplitudes involved in Eq.(71), so 
that \T semit t t j\ < dim (H) = N is always satisfied. These cancellations are due to the 

leading complex terms exp ^— i^j^j- 

• In some specific examples, namely the linear map M , Eq.(2), with particular values 
of 27r/l, we observed in [27] that at t ~ 2t El all the actions exp ^— i 1 ^- j = 1 are equal 

and add together. This implies that the upper bound Tr (j&taruf) = dim (H) = N 

is reached, and therefore that M^rus 00 Id which implies revivals of quantum states 
and existence of scarred eigenstates, i.e. non equidistributed invariant semiclassical 
measures. 

• Let us give now a "philosophical" remark which shows again that t = 2t E seems to 
be a critical value of time. If we define the "complexity" of the trace formula to be 
equal to the number of periodic orbits, it grows like e Xot . This complexity is larger 
than the "complexity" of the linear quantum problem (equal to the size of the matrix 
M: N 2 ~ e 2A ° iB ), for t > 2t E . Again this shows that for t > 2t E , periodic orbits 
manifest themselves in the semiclassical trace formula through collective averaging 
effects (see Berry in [29], or [25] where this open problem is also discussed). 

6.2 Numerical results and observation of self-averaging effects 

We illustrate the semi-classical formula for the trace Eq.(70), with the example defined by 
Eq.(12). We have computed Tr ^M*j numerically for t = — > 11 (this is easy because 
the Hilbert space is finite dimensional, and we expect the numerical result to have a good 
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accuracy). We compare it with the semi-classical approximation T se mi,t,j, Eq.(71), where 
every normal forms has been computed numerically 6 up to a given order J > 2. We choose 
successively N = 1/ (2nh) = 10, TV = 1/ (2nH) = 100 with J = 2 and J = 4. The 



observed error is denoted by Errorj (t) 



dcf 



The practical limitation 



Tr yM 1 — T semi tt j 

in this numerical calculation is the increasing number Nt — e Ao< of periodic orbits; already 
Mt=n = 39601. In Theorem 12, we have estimated that 



with ej (t) = h K = h^ c = n J ' 2 - t ' tE 

with tfi 

On Figure 7, we plot three functions of (i/is) on a logarithmic scale: 



Errorj (t) < C X Ej (t) 
tlog(lM). 



Tr M 



the 



observed error Errorj (t) and the upper bound of the error ej (t). We also draw the line 



Tr 



Max 



1/ (2irh) (because obviously 



< dim ft 



torus 



1/ (2%h) ), and the line of 



Tr [M 1 

the particular value e%j = h J l 2 which is expected to be an upper bound of the error from 
heuristic arguments given below. The unexpected observed fact is that Errorj (t) < e^j, 
although Errorj (t) < Ej (t) = Sf iy .jh~ t ^ E has only been proved. 



Remarks and observations: 



M is an unitary operator, so 



Tr M 



< dimHtoms = 1/ (2nh). However our error 

estimate gives ej (t) = h J ^ 2 ~ t ^ tE > 1/ (27rh) for t > + l)t E . So our estimates is 
not interesting fort> + 1) t E - 



According to Eq.(69), the error term Errorj (t) is a sum of Af t 



3 Aot 



complex num- 



bers, each being bounded in magnitude by B ~ e - x ot/2^j/2 ^^jg appeared in the 
proof of Proposition 11). With an heuristic point of view, if we think that these num- 
bers behave as independent random variables, then the total error is estimated (by 
the central limit theorem) to be of order \fW t B ~ h J ^ 2 as the numerical calculations 
shows. It would a nice result to explain such a behaviour. 



7 Conclusion 

The main result of this paper is the use of semiclassical non stationary normal form de- 
scription of the hyperbolic dynamics first introduced by David DeLatte [19], and which 
gives some invariant semi-classical cocycles of the hyperbolic dynamics. With them, we 
obtain semi-classical expressions for long time of order t ~ Clog (l/h), with any C > 0. 
These cocycles are a generalisation to all order in non linearities and in power of h of the 
Lyapounov cocycle which controls the linear instability of the dynamics. These cocycles 
allow us to have a control on the error term of the trace formula (or other semi-classical 

6 The algorithm to compute the normal forms is explained in Section E. 
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ft = 0.1, J = 2 




ft = 0.1, J = 4 




0.2 0.4 0.6 0.3 1 1.2 1.4 1.6 1.3 2 2. 

ft = 0.01, J = 2 



tftE 




0.2 0.4 0.6 0.3 1 1J 1.4 1.6 1.3 2 2. 

ft = 0.01, J = 4 



Figure 7: Numerical results for 



Tr [ M* 



and its semiclassical approximation T semi ^ j. 



We have plotted Errorj (t) 



clef 



Tr M 



semi,t,J 



in log. scale. 
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formula) for long time. We could ask the question if this Non-stationary Normal descrip- 
tion presented in Section 5 is really necessary? the author believes that it is, because it 
gives Eq.(45) where T x depends on the point x and not on the trajectory x (t) = M t x. 
Without this result, we could expect that the Normal Form construction would depend 
on the trajectory x (t) as a whole, and could diverge with t in an uncontrollable way. We 
would be then blocked. 

There are many perspectives suggested by this work. Numerical calculations presented 
in Section 6.2 seem to show that our bounds are not sharp, and that the actual errors 
are much smaller. In particular, Gutzwiller trace formula at the leading order J = 2, i.e. 
without semi-classical normal forms corrections, could be correct for these "long time", with 
an error O (K) uniform in time. To understand this surprising fact we should understand the 
way all the complex semi-classical contributions add and compensate together. We think 
that the Ruelle-Pollicott thermodynamical formalism of transfer operators in the context 
of prequantum dynamics, could help us in that direction. Some preliminary results in this 
direction are obtained in [25], in the case of a linear hyperbolic map, and we hope to be 
able to extend these results to non linear maps. 

An other perspective is a generalization of this approach to a uniform hyperbolic flow, 
like the geodesic flow on a negative curvature manifold. We think that the approach of 
"semi-classical non stationary normal forms" developed in appendix D could be generalized 
for such models. A similar but less obvious generalization could be done for maximal 
hyperbolic sets, situations which are met in general mixed chaotic systems, with homoclinic 
intersections of stable/unstable manifolds giving horseshoes [45]. 

A Proof of Theorem 6 page 22 on semi-classical evolu- 
tion in a neighborhood of a classical trajectory for 
long time 

We first introduce some notations specific to this appendix. 

A.l Some notations 

Let D C M 2 (a measurable set). We define: 




where \x) is a coherent state defined by Eq. (25). One has the closure relation: 



Pd + Ps*\d = Id 



(72) 



Proof. Id = /j 



R 2 2nh 



x)(x\ = /j 



D 2-itH 



+Ji 



R 2 \D 2-irh 



x) (x 



□ 
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Let < (5 < 1/2 (The sharper results will correspond to j3 close to 0). Let D C M. 2 a 
measurable set, or a sequence of sets D^ which depend on h, with h —> 0. In particular .D 
may be a point D = {x} } x G M 2 . 



Definition 13. Let's define the "exterior domain" E(D), by: 



E (D) = {x G M 2 such that dist (x, £>) > fr 1 / 2 ^} 



(73) 



and its complementary, the "interior" I (D) 



1(D) = R 2 \E (D) 



Remarks: 1(D) is just the domain D which has been thickened by h 1 ^ 2 13 . In particular 
D C 1(D). In Eq.(39), we have defined P x ,/3, which in the present notation is Pi{ x y. 



In this appendix, the classical dynamics is the map M : M 2 
and the quantum map M has been defined in Eq.(15). 



I 2 defined in Section 2, 



A. 2 Semi-classical evolution after finite time 

The results of this appendix will rely on the following lemma which is a quite standard 
result on semi-classical evolution of wave packets after finite time (with respect to h — * 0). 

Lemma 2. For any K G N, and H small enough, there exists C > 0, such that for any 
x G M 2 , and any domain DcK 2 such that D C E (Mi), 



P D M\x) 



L- 



< Ch 



K 



(74) 



where \x) is a coherent state defined in Eg. (25). 



Remarks 



In all this appendix, we will use this result, and write: e = Ch . At the end, K will be 
chosen large enough. In common notations, we can write that PdM\x) = O (h°°) 

L 2 

uniformly over x G M 2 and D C E (Mx) . Lemma 2 means that the evolved state 
M\x) is "localized" at point Mx at order O (h°°). In our case, the uniformity in x is 
due to the fact that the dynamics is defined on a torus (compact space). 

Lemma 2 follows directly from more precise results obtained by A. Joye and G. 
Hagedorn in ref. [36] (theorem 3.2), or M. Combescure and D. Robert [13], theorem 
3.1, or A. Iantchenko [41], Lemma 5. 



Corollary 14. Suppose that the domain D C R 2 has finite measure S(D) = f f x£D dx. 
Then for any domain D' C E (M (D)) 



Pd ,mp d <mf\ 

L 2 \ 2ixh J 

(We don't think that the dependence on D of the right hand side is sharp). 
Proof. Let ip G L 2 (R) normalized, then 

d 2 x 



(75) 



P D ,MP D ^ 



L 2 



< 



< e 



P D ,M 
d 2 x 



2ttH 



x) (x\i/j) 



L 2 



D 2nh 
d 2 x 



|(a#)| P D >M\x) 



2irh 



IW)l<(^V /2 e 



In the last line we have used Eq.(74) because D' C E (M (I?)) C E (Mx) if x e D. We 
have also used Cauchy-Schwarz inequality. □ 

A. 3 Long time Semi-classical evolution 

Let D C M. 2 be a set at time t = 0, with finite measure. Let 

h (D) D 

And for any t > 1, 

3t (£>) = £ (M/t-O , I* (£>) = M 2 \^ (D) 
I t (D) is just obtained from .D, by evolution and thickening at each step, see figure 8. 

E 2 {D) 




etc... 



Figure 8: Description of the sets E t (D) and It(D), constructed inductively from D, by 
evolution and thickening. 
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Lemma 3. For any time t > 1: 



2Tih J 



L 2 



< e 



,s=0 



l/2> 



(76) 



Proo/. From Eq.(72), we get Id = P Et[D) + P It(D) , and E t (D) = E x {I t _ x (D)). Eq.(75) 
gives for any t: 



? El{D) MP D 



< 



l 2 \ 2nh 

where e does not depend on t G f , and D. By induction on time t we get: 



P Et{D) M t P D 



L 2 



P Et(D) M (P Et _ l(D) + P h _ l{D) M- L P D 



L 2 



< \\P Et (D)M [P Et _ liD) M^P D 



L 2 



+ 



P EtiD) MP It _ l(D) ) M^Po 



L 2 



1/2 



where we have used: |v4 + £>| i2 < \Ml 2 + I-^Il 2 ' I^-^Il 2 — I^Il 2 I-^Il 2 ' an d 



P 



D 



M 



L 2 



L 2 



1. 



< 1, 



□ 



A. 3.1 A particular application 

We consider now a particular application of the previous results, which will be useful for 
hyperbolic dynamics. Let t > 1, and suppose that Do, A, . . . , A is a sequence of sets. For 
any s G [1, t — 1], let 

G s = h (A-i) \A 
which corresponds to points coming from A-i but not belonging to D s . 

Lemma 4. Suppose that for every s G — 1], 

J t _ s (G s ) n A = 

which means that "points which leave the domain D s do not come back to the domain A " 
(despite the thickening procedure), see Figure 9. Then 



Pd.M'Pdo - P Dt MP Dt ^ ■ ■ ■ MP Dl MP Do 



L 2 



t-l 



< 



2ttH 



/2 



fc=0 
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Proof. For any s 6 [1, t — 1], let 

F s = E 1 (£>._!) \£>„ G s = /! (£)._!) \D S 

So that I 2 = D s U F s U G s as a disjoint union. The hypothesis I t - S (G s ) H D t = $ means 
that D t C E t _ s (G s ). Then Eq.(76) gives, 



< e 



't-s-l 

E 

fc=0 



S(4(G S )) 



2nh 



l/2> 



Also, P s C -Ei (As-i), so Eq.(75) gives: 



< 



L 2 



S(D a ) 
2nh 



1/2 



Now we proceed by induction on s. We use Id = Pd 1 + Pf\ + Pgi (from the disjoint union 
M 2 = D, U F, U G,V and write 



PB t M l P D , = Pd.M*' 1 ( P Dl + P Fl + P Gl ) MP Do 



P Dt M*P Do - PdM^PuMP, 



Dl ivir Do 



< 



L 2 



Pd.M 1 ' 1 [P Fl MP Do 



L 2 



PdM^P^ 1 MP; 



Do 



L 2 



< e 



1/2 



i-2 

E 

,fe=0 



S(/ fc (Gi)) 

27Tk 



1/2' 
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Similarly, 



P^M^PnMPn, - PnM^PnMPnMP, 



D 2 ivi r Dl ivi r Do 



L 2 < £ \ 

't-3 



S(L> 2 )^ 1/2 



E 

,fc=0 



S(4(G 2 ))\ 1/2 



2irh 



and therefore after few similar steps 
Pd.M'Pd, - P Dt MP Dt _ x . . . P Dl MP Do 



L- 



< e 



t-i 

E 



2ttH 



1/2 



't-1 t-s-1 

EE 

,s=l k=0 



S{h(G a )) 
2nh 



l/2> 



□ 



A. 4 Application to hyperbolic dynamics 
A. 4.1 Proof of theorem 6 

The preceding results were rather general. We will apply them to the hyperbolic dynamics 
considered in this paper. In order to prove theorem 6, we consider a trajectory x (s) = 
M s (x (0)), s = 1 . . . t, and disks D s of radius ft 1 / 2 "" at each point x (s). Using definitions 
introduced with Lemma 4 page 40, we remark that the hypothesis I t - S (G s ) D D t — is 
satisfied, because if j3 < a, and due to hyperbolic instability of the trajectory x(s), the 
dynamics sends the domain G s away from the domain D t , in spite of the thickening (see 
figure 9). We have S (D s ) ~ O (fi 1 / 2-0 ) but the surface S(Ik(G s )) grows exponentially 
fast with k (in the unstable direction) because of the thickening at each time step. A rough 
estimate gives: 

S (4 (G)) = O (( e Xma * k h 1/2 - a ) (h 1/2 ~ a )) 

Suppose that < t < C\og(l/H) e Xmaxt < frr CXmax . Then Lemma 4, and e = Ch K , 
gives 



P Dt M l PD - PD t MP Dt ^ ■ ■ ■ MP Dl MP Do 

K t-l / t-s-l 



L 2 



Ch 



< 



C 



k=0 

h K h- a h- CXmax/2 t 



For any iiT' > 1, the bound is O (fo K ) if K is chosen large enough. So we get Eq.(41). 
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A. 4. 2 An other consequence 

Lemma 5. Let < a < 1/2. For any C > 0. any t, such that \t\ < C\og(l/h), for any 
n E 1? , any x E M 2 , such that \x — x n j\ > h 1,/2 ~ a (which means that x is at some distance 
of the hyperbolic fixed point x n ± of the hyperbolic map M n ^, Eq.(34)) then 

(x\M n , t \x) = 0(h°°) (77) 

(uniformly other x E M. 2 ,n,t), where M n j = T- n M l has been defined in Eq.(33). A light 
generalization is 

(x'\M n<t \x) = O (ft 00 ) (78) 
if \x — x Ht t\ > K x l 2 ~ a and \x — x'\ < h 1 / 2 ~ a ' , with < a' < a. This implies: 

Tr (P Xntt , a M n , t P w ) = Tr [M n , t P Xn ^ + O (ft 00 ) (79) 
where P Xn , t ,a is defined in Eq.(39). 

Proof. Let < 7 < a', and let D' x be the disk of radius ft 1 / 2-7 with center x. With the 
hypothesis \x — x nit \ > ft 1 / 2- ", and \x — x'\ < H l ^ 2 ~ a , one checks that for < (3 < 7, then 
D' x , C E t (D' x ). Then from Lemma 3, and similar estimates as in Section A. 4.1, one gets 



P D >M n +P D , 



0(h°°). This implies that (x'\M n>t \x) = 0(h°°). To show Eq.(79), write 



Tr (Mn.tP^.a) = Tr (P w M n , t P Xn ^ + Tr ((l - P Xn ^ M n , t P Xn ^. But 

jHf Tr (h _ p \ M n>t p\ = t [ ^- 2 { X \x'){x'\M n , t \x) 



< 1, 



Using the exponential decrease of = C exp (— \x — x'\ /fx), and (x'\M njt \x) 

one first obtains: 

T= [ [ ^^(x\x')(x'\M nit \x) + O (h°°) 

with < o>i,a!' < a, and using Eq.(78), one finally obtains X = 0(h°°), and deduces 
Eq.(79). □ 

B Proof of Lemma 1 page 23 for the trace of an indi- 
vidual orbit for long time 

B.l General results to control localization of wave packets 

We first give some results which will allow us to control the escape towards infinity in 
R 2 of the phase-space distribution of wave packets after long time, under a hyperbolic 
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dynamics. The usual semiclassical Egorov Theorem which describes the propagation of 
quantum observables does not work after a fixed constant times the Ehrenfest time (He 
in our case, see Section 1.2). The aim of the Section is to develop a weaker description 
which will be able to reach t = C\og(l/h), with any C > 0. Remind that the classical 
phase space is x = (q,p) G IR 2 . The quantum Hilbert space is Tipiane = L 2 (R) with Planck 
constant h. 

Let m be an order function which possibly decreases at infinity ([21] p. 81, [49] p. 12, 
[23] p. 52). 



Definition 15. A weight function is a semi-classical symbol W : M 2 
S (m), such that = max^ (W (x)) = 1, and W is elliptic. 

Remark: W Elliptic means \W\ > jjm, with C > 0, then W^ 1 G S 
p.100, [49] p. 13. 



m 



, W G 
- 1 ), see [21] 



Example: W (x) = m(x) = 1/ (x) with (x) = y/l + q 2 + p 2 , k > 1, or W (x) = 
m (x) = exp (— (x)). 

We write W = Op Weyl (W). We check that these operators have dense domain on 
L 2 (R). 



Definition 16. If W is a weight function, a sequence of states iph, with h 
localized if ||^|| = 1 and if there exists C > (independent of h) such that 



0, is W- 



L 2 



< c 



The intuitive idea of this definition, is that a W-localized quantum state has a distribu- 
tion in phase space which is bounded by the function W (x). This definition is very similar 
with weight functions used to treat tunnelling effect. 

Remind that the dynamics considered in this paper is the map M : M, 2 — > M 2 , defined in 
Eq.(3) as a product M = MiM , with M G SL(2,R) linear hyperbolic and M x bounded. 
We also defined M in Eq.(15). 

Proposition 17. Suppose that W G S(m ), W\ G S (mi) are weight functions such that 

Wx (x) > W (M~ l (x)) , Vx G M 2 
and mi (x) > m (Mq 1 (x)) , \/x. Suppose that ipo is Wo-localized. Then 

fa = M^o 

is W\-localized. 
Proof. Write 

From Egorov theorem ([21] p. 125., [49] page 138, [23] p. 139) 

MWoM- 1 = w' + hk 
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where W' = Op(W) and W has Weyl symbol W (x) = W (M~ x (x)), and belongs to 
S(m'), with order function m! = moM^ 1 . The remainder R has Weyl symbol R G S(m'). 
Then 



W^MxWWo'Vo = W{ 1 (W' + fri?) MWq-Vi 

Wf 1 (#* + ft£ 



W -Vo 

By hypothesis, Wi > W & l^f^'L < 1, and Wf 1 (W + hR) e S (m^rri) = S(l). 
By Calderon-Vaillancourt Theorem ([21] p.85, [49] page 43, [23] p.69), W^ 1 (w' + hR 



< 



l + 0(h). Also 



M 



1, and 



< C. We get: 



WTVi <c.(i + £>(£)) 



□ 



Proposition 18. Suppose that Wq , W\ are weight functions, and that ipo is W^-localized, 
and ipi is W\-localized. Then there exists C > and h > such that for any H < Hq: 



Proof. We suppose 



< Co and 



< C x . Then 



< 



From Calderon-Vaillancourt Theorem, 



W W 1 



< {WoW^ + Oih). We get: 



But |WWiloo < 1, so 
with C > 0. 



IMM>| ^CodiwWiL + o^) 
l(^i>l <c*.|wWiL 



□ 



B.2 Application 

The previous paragraph will allow us to obtain the following Lemma: 

Lemma 6. There exists R > 0, such that for any C > 0, < a < 1/2, there exists Hq, C\, 
such that for any t < C/h, h < h , \xq\ > R, one has 

(M--R) 

(x |M*|x ) < Cie 
where \xq) is a coherent state at position Xq = (qo,Po) £ K 2 - 
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This last results means that the auto-correlation function 
nentially fast as xq goes to infinity. This is due to the hyper 
M on R 2 . The proof is given below. We deduce: 



(xqIM^o) 



decreases expo- 



Dolicity of the classical map 



Corollary 19. There exists R > 0, C > 0, such that for any 1 < t < C/h, K < h Q , one 
has 



(x |M*|z ) 



|x |>.R 



(x\M n;t \x) 



dxo 
2ttK ' 

dx 



\x — X n ,t\>R 



2nh 



o(h°°] 

= 0{ff 



(80) 



The second equation is a easy generalisation of the first one, where we considered the 
map Eq. (33) . This last result improves the result of convergence Eq. (36) , because it shows 
that the integral is semiclassically negligible outside a disk of fixed radius R. 



B.2.1 Proof of Lemma 6 

In all the paragraph, let us denote by x = (q,p) £ M 2 the unstable/stable coordinates in 
which the matrix M , Eq.(2), is diagonal. M : (q,p) — > (e x °q, e~ x °p) . 

A first result on the classical map M: Remind that the perturbation Mi is bounded: 
| Mi (x) — x\ < C, for any x £ K 2 . Without further assumption on Mi (we don't assume 
here that M = M\Mq is hyperbolic), we have the following simple result Lemma on the 
map M = MiM : 

Lemma 7. For any e, such that < e < A , let X = A — e > 0, R = ff e -°) > 0. For any 
point x = (q,p), let x' = (q',p') — M (x). If q > R, then q' > e x q. 

This result guaranties that under the map M, points escape exponentially fast towards 
infinity if q is large enough. 

Proof We have x 1 = M 1 M (x), so \q' - e x °q\ < C, so q' > e x °q - C. And e Xo q - C > 
e x q^ q > R. □ 



Choice of Weight functions W t : In order to apply later on Proposition 17, let x = 
(QoiPo) ^ with q > R. Let < a < 1/2 (by rescaling, we use semiclassical calculus 
at the scale h l / 2 ~ a 1 [21] p. 82, [23] p. 52). For any t £ N, let us choose a weight function 
W t (q) (independent of p and increasing in q), defined by W t (q) = 1, for q > qoe xt , W t (q) = 

exp ( (qe~ xt - q ) /h 1 / 2 ^) for e xt R < q < q e xt , and W t (q) = W.^ = e ^ q °- R ^ hl/2 ~ a for 
q < e xt R. We smooth W t (q) in the vicinity of R and q and let it be still increasing but 
C°°. For any t, the function W t is a weight function according to definition 15, but at the 
semiclassical scale x = xj (Ti 1 / 2- ") . From Lemma 7, we check that it satisfies: 

Lemma 8. For any x = (q,p) £ M 2 , any t > 1, then 
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W t (x) > W t -i (M" 1 (z)) 

Lemma 9. Let ipo — \ x o) be a coherent state at position x G M. 2 , with go > R. Then 
tpo is W^-localized (according to definition 16), and for any < t < C/h, ifj t = f M^o « s 
Wt-localized. 

Proof. ipQ is Wo-localized because as explained in Section 2.2.3, a Gaussian wave packet 
ipo (g) is exponentially localized in q, around go- Then with Lemma 8 and Lemma 17, we 
deduce iteratively that iftt is ^-localized for any t. We have to be careful, and observe that 



in the proof of Lemma 17, there appears 



< C t -i- (1 + (K)). Under iterations, 



this gives C (l + O(h))\ which is still O (1) if t < C/h. □ 

Last step of the proof: Now let us consider another weight function W (q) (again 
independent of p) defined by W (q) = 1 for q < q , W (q) = exp (— (q — q ) /ft 1 / 2 "") if 
q > q . We modify W in the vicinity of go, so that it is C°°. As above, it is clear that 
the coherent state tpo = \xo) is W-localized. Then we apply Proposition 18 to deduce that 
KV»<M)| < C. IWWtl^. Observe that {WW^ = W.^ = e -(io-R)/h^ so this giyes 

(xolM'lxo) < Ce^ R)/hl/2 ~ a (81) 

The same analysis can be done for g < R, and similarly for \p \ > R (in that last case, 
this is the stable direction, we have to work in the past, writing (x |M*|x ) = (^_t|^>o))- 
Finally we deduce Lemma 6. 

B.3 A more refined estimate 

Corollary 19 is not precise enough to give Lemma 1 we are looking for. Indeed, instead of 
a disk of finite radius R, we have to reduce the integral Eq.(36) to a disk of smaller radius 
h 1,2 ~ a . Let < a < 1/2, and t < Clog(l/ft), with C > 0. Decompose Eq.(36) in three 
parts: 

T(M„, t ) = / (x\M n , t \x)- d: 

V / J\x-Xr, t\<h 1 / 2 - a 



2ivh 



f •*> dx f a dx 

+ / (x\M n , t \x)—r+ / (x|M n)t |x)-— 

Jh 1 / 2 -°'<\x-Xn,t\<R n JR<\x-x n ,t\ Z7ltl 

From Eq.(80) and Eq.(77), the last two integrals are 0{h°°). One obtains 
T(M n , t ) = [ ^-{x\M n>t \x) + O (h°°) = Tr (M nt P Xnt ^\ + O (h°°) 

where P Xn>t is defined in Eq.(39). Eq.(79) gives Tr (m^P^,,) = Tr (P Xrht M n ,tPx n , t ) + 
O (h°°). Finally Lemma 1 is proved. 
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B.3.1 Estimates of Trace norm in terms of Operator norm 

We consider operators P XtCt denned in Eq.(39), with x — 0. 



Lemma 10. Suppose that A is bounded: 
K > 0. Let a' > a. Then in Trace norm: 



A 



L 2 



0(1), and 



P a A 



L 2 



o(h K ) with 



Similarly if 



D 



P n >A 



L 2 



Oh 



-K-a' 



P a /AP a i 



Oh 



-K-a' 



0(1), and 



P a BP a 



L 2 



0{h K ) then 



P a tBP a r 

Proof. We first give basic estimates: 



Oh 



-K-a! 







Pa 





If a' > a, 



and 



One has 



\x\<ro-/< 



P n I 1 — P n i 



-a 2llh 



(h- 2a ) 



(82) 



(83) 



L 2 



o{h°°) 



Pn 



L 2 



< 1 



Pa' PaA 



< 



1 = 1 - P a + P a , and 



where 

Pa> (l ~ Pa ) 

Therefore P n >A 



Pa' 



Pa'A 



PaA l 



L 2 



(h K ~ a ') (see [30], prop 5.4 page 62). Write 



< 



Pa' PaA 



Pa' ( 1 Pa I A 



< 



Pa' 1 - Pc 



A 



< 



L 2 



Pa' [I- Pc 



L 2 



A 



L 2 



0{W 



0(h K ~ a ') and 



Pa'AP a i 



the second estimate, let A = BP a - This gives 



< 



Pa'A 



Pa'BP a 



Pa' 

L 2 

Pa'A 



1 



above, one can show that 



P a ' BP a ' 



< 



P a 'BP a 



o{h c 



0(h K ~ a '). For 
: (h K ~ a '). As 
□ 



C Useful results for finite time evolution 

In this appendix we collect well known results on semiclassical finite time evolution. We 
state them in the context of the present paper. They are used in other parts of the paper. 
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C.l Truncation of Taylor series of the Hamiltonian 



Theorem 20. Let Hi (x,t) and H 2 (x,t) be two symbols on x = (q,p) G M?, such that 
is a fixed point (i.e. (DHi) x=0 = 0), and suppose that they have identical Taylor series in 
(x, h) at the origin (0, 0, 0), up to order J: 



j+i 



(i.e. Taylor terms fr!q a p b are identical if I + a + b < J). Let U\ = exp (—iHi/h\ and 

U 2 = exp (^—iH 2 /hj , with Hj = Opweyi(Hj), j = 1,2. Let < a < 1/2, and P a be 
defined by Eq. (39). Then 

< Ch A (84) 



with A — (I — a') ( J + 1) — 1, for any a' > a. 



This result is interesting in the semiclassical limit if A > <^ J > Yzff ? f° r example if 
J = 2 and < a < 1/6. This theorem will be used page 61 for the proof of Theorem 7. 

C.2 Proof of Theorem 20 

We will use the Duhamel formula [44] (which is not semiclassical) : 

Proposition 21. Suppose that H (t) is a self-adjoint operator in a Hilbert space 7i, and 
"0 (t) 6 H is solution of 

where £ (t) G H. ||£ (t)\\ < \i (t), and \\ip (0)|| = 1. Suppose that ip (t) is solution of 

.dip 



with cp (0) = if) (0). Then 



U{t)-ip{U)\\<\f\{s)ds 



(85) 



We can deduce a "semiclassical version" of the Duhamel formula: 

Corollary 22. IfHi,H 2 G S'(l) are real bounded semiclassical symbols (cf [21], p. 81), 
and if K = H l - H 2 G S~ k (1) (i.e. K G h k S{\)) , and U x = exp (-iHjhY U 2 = 

exp ( —iH 2 /h), with Hj = Opweyi (Hj), j = 1,2, then 



U x -U 2 



L' 2 
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A 



L 2 



Proof. A e S- k (l), thus 
if) G H, be normalized: 

with ipi (0) = if) 2 (0) = if; . This gives ihdfa/dt 



O (h k ) , from Calderon-Vaillancourt Theorem. Let 



1, and if>i (t), ip 2 (t), solutions of ih^- = Hjif>j, j = 1,2, 

#1 - A) ^2 = #1^2 + £ (t), With 



!l£(*)ll< 



A' 



L 2 



C (ft fc ). From Eq.(85), we deduce 



U 2 if) - 



L 2 



<\0{h k ). □ 



Micro-Localization Lemma: 

Lemma 11. Let < a < a' < 1/2, and H (x,t) a real symbol, with fixed point 0: 
DH/ X=Q = 0. Let H (t) = Opweyi (H (t)), and H' (t) = H (t) P a > be the truncated Hamilto- 
nian, where P a >is defined in Eq.(39). Let U (t) and U' (i) be the evolution operators defined 

by = H (t)U {t), ih^^ = H' (t) U' (t), and U (0) = U' (0) = Id. Then for any 

finite time t, 



U{t)-U' (t))P a 



L 2 



o(h c 



Remarks: Remind that P a truncates on a disk of radius ft 1 / 2 " , and notice that ft 1 / 2 <C 
^ f^l 2 ~ a in the semiclassical limit. Proposition 11 means that the evolution of a 
quantum states in the vicinity of a fixed point (or more generally of a trajectory) after finite 
time, depends only of H (x) in the vicinity of this fixed point. The proof of Proposition 11 
relies on Egorov Theorem, and can be found for example in [41], Lemma 3, Lemma 4. 



Proof of Theorem 20: 

Proof. We use notations of Lemma 11 and Corollary 22, and combine them. Let H'- = 
HjX'i j = 1)2, where %' (x) truncates at radius r' = h l ^ 2 ~ a ' (i.e. %' (x) = 1, if |x| < r', 
x' (x) = 0, if \x\ > 2r'). Similarly, let x i x ) which truncates at radius r = h l / 2 ~ a . Let Uj = 
exp t-iHj/h), U'j = exp (-ify/h). By hypothesis, we have A = (H 2 — H[) e S~, k (1), 



with h k = h^ 2 -^ i.e. k = (1/2 - a') J. From Corollary 22, U[ - U' 2 
Let x = Op Wey i (x). One has 



L 2 



o{h k - 1 ) l 



U 1 -U 2 )x 



L 2 



- U[ + U[ - U' 2 



u> 2 



U2)X 



< 



Ui-u[)x + (u'l-utix 



U'o-U, 



X 



< o + o (h k ~ 1 ) + o (h c 



□ 
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D Proof of Theorem 7 on Semiclassical Non Stationary 
Normal Forms 

In this appendix we develop the Theory of Semi-classical Non-Stationary Normal form for 
hyperbolic map on the torus, which extends the results of David DeLatte [19] from the 
classical to the semi-classical case. 

D.l Description of the dynamics: 

Let us write x = (q,p) G R 2 . Remind that the classical map M = MiM : R 2 — > R 2 was 
defined in Eq.(3) as the result of a first flow M generated by the quadratic Hamiltonian H 
for time interval t G [0, 1], followed by the "perturbation" Mi generated by the Hamiltonian 
Hi on time interval t G [1,2]. Hi (x) is periodic with respect to Z 2 C R 2 . So the dynamics 
is defined by a time dependant Hamiltonian function written H (x, t) with period 2 in t. 

For any t, t' G R, we denote M t>t i : R 2 — >• R 2 the flow generated by the Hamiltonian 
H(x,t), in the time interval [t, t']. In particular M = M 0j 2- 

Using the Weyl quantization procedure, we have defined the corresponding quantum 
Hamiltonian H (t) = Opweyi (H (x,t)), in Section 2.2. Conversely H (x,t) is called the 
total semiclassical symbol of H{t). The unitary evolution operator M t ,t? is defined by 

^ = -iH (f) M t ,t', and M t ,t = Id. 

Stable and unstable tangent vectors for intermediate time t: In Section 2.1.3, we 
explained that for every point x G T 2 , there is a basis of tangent vectors (u x , s x ) G T X T 2 , 
continuous with respect to x G T 2 , and tangent to the unstable/stable foliation. We can 
define this basis for intermediate time t G [0, 2] (and then for any time t G R by requiring 
periodicity), giving vectors u Xt t,s x j as follows: the direction of u Xjt (respect. s Xt t) is the 
direction (DM 0tt ) x , u x > with x = M 0jt (a/) (respect, for s x r), and we fix the choice of u x>t 
(respect. s x j) by requiring that 

u X)t A s X)t = 1 (86) 

i.e. they form a symplectic basis, and we also impose that Hwas.tH — ll s x,t||- 

Let us define the lattices r = Z 2 , r t = M 0>t (T ) for t G [0, 1], and T t = T for t G [1, 2]. 

Let us define T t = f M. 2 /T t . If we consider i G S 1 = f / (2Z), then the phase space and 
time containing points (x, t) forms a non trivial bundle of tori B — > S 1 over S , whose fiber 
over t G S 1 is the torus T 2 = M 2 /r t . u X)t and s Xjt are continuous in x G M 2 and periodic 
with respect to T t , and also continuous and periodic in t G R. In other words they are 
continuous functions of (x, t) G B. In this section periodicity in x and i will always refer 
to the lattice T t and to S 1 = R/ (2Z) respectively. 
Define the expansion rate A (x, t, t') by 

(DM t , v ) x (u x , t ) = e A ^'W t , ( ,(*),t' 
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and the local expansion rate rj (x, t) by 

r) (x, t) is a piece-wise continuous function and periodic. In particular 

X x —[ V (x (t) , t) dt = X (x, 0, 2) 
Jo 

gives the expansion rate already defined in Eq.(7). 

D.2 Preliminaries on Time-dependant canonical transformations 

In order to simplify the description of the dynamics, and obtain the non-stationary normal 
form, we will perform time dependant canonical transformations at the semi-classical level 
([49], [23], [21]). In this paragraph, we prepare some results and fix the notations. 

D.2.1 Moyal formula 

We remind that if A, B are two symbols and A = Opweyi {A) and B = Opw ey i {B), then 
the Moyal product Jj between symbols correspond to product of operators and is defined 
by ([49] p. 41) 

A$B = a Weyl (ABj (88) 

Where aweyi gives the symbol of an operator. The Moyal formula gives the symbol of a 
product of operators: 

AlB = m n ^7 (AV n B) = AB + ih\ {A, £?} + ... (89) 



n>0 



where V = 4-4- - #-#-, and {A, B} = (AVB) = §^ - §^ is the Poisson Bracket. 

op aq aq ap 1 L 5 j v j op oq oq op 

As a consequence the symbol of the commutator is given by 

^ (B) = *w eyl Q [A, B])=J2 (iHr 22n+1 " - 1V (AV^B) (90) 

\ ' n>0 ^ 

= {A,B} + ^(AV*B)+... 
The first term, which is the classical Poisson Bracket will also be written: 

%class,A (B) d = {A, B} 

These formula exist in a more general framework of star-product, and all what follows 
can be applied to any star-product, in particular to any choice of quantization procedure. 
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D.2.2 Time dependant semiclassical canonical transformation 

Let G t (x) be a semiclassical symbol which depends on t, G t = Op Weyl (G t (x)) and U t = 

exp y—iGt/nj the unitary transformation in L? (R) generated by G t for any fixed value of 
t. 



Proposition 23. For any (t,t'), the transformed evolution operator: 

M' t>t , = U^MtjUt (91) 

is generated by a transformed Hamiltonian symbol H' . (i.e. dM' tt ,/dt' = — ^H'{t')M' tt , } 
and M' t t = Id, with H' = Op Weyl (H 1 ) ), given by 



H' = exp(-X Gt )H + R (92) 



with: 



o \- yfc fdGA fdGA 1 fdG 



(fc + 1)! Ut v 5t ; V dt J 2 " V dt 



Notice that the term R is due to the time dependence of the canonical transformation, 
and that exp (— X c iass,G) H — H o U^ 1 is the function H after the classical canonical 
transformation generated by G t . 

Proof. For the proof, we consider the enlarged phase space (t,r,q,p) G R 4 with the sym- 
plectic two form ui = dq Adp + dt A dr, where t is considered here as a dynamical variable, 
and r is its conjugated variable. This is a well known and useful trick both in classical and 
quantum mechanics which allows to map a time dependant dynamics onto a time indepen- 
dent dynamics, so that we can use well known tools. With f = —ihd/dt, and H = r + H, 
the time-dependant Schrodinger equation ifr^ = Hip can be written 

H^> = (f + H\vj) = (93) 

i.e. as the time-independent Schrodinger equation for ip(q,t). We will use bold fonts 
when dealing with the enlarged phase space or enlarged Hilbert space L? (R 2 ). Let G = 

Opweyi (G), where G (q,p,t,r) = G t {q,p) is independent of r, let U = exp ^— iG/fij, and 
consider the transformed Hamiltonian 

H' = U^HU (94) 

whose total symbol is: 

H' = exp (-X G ) H = H - X G (H) + ^X 2 G (H) + . . . 
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where X G (H) = \ (e,> (-1)* j^yG (f f +1 H ) , with D = (J | - g *) + (|| - ||) = 
I?i + P g . Because G does not depend on r, we have GDH = (— ^) + GT> g H = 
(-ff ) + G t P ? # and for k > 2, G D fc H = G (V t + £> g ) fc H = G t V k q H. Then X G (H) = 
(-ff ) + £ Gt (#), and X k G (H) = -S^; 1 (^) + X k Gt (H), for fc > 1. Therefore we get 

H' = exp (-X G ) H — t + H' 

with 

#' = exp(-X Gt )H + R 

Now, for any u G R, let M u = exp ^-mH/sj and = exp (— iuil'/?i\. From 
Eq.(94), we have the conjugation relation M' u = U^M^U. But (q, t) = (& t p \ (q, t) 

and (wiuip^ (q,t) = \Mt-u,ttyj (q,t — u )- So the conjugation relation can be written: 
MU >t = U^Mt-uA-u, or M[ t , = U^MtjUt, for any t, t' . □ 



D.2.3 Translations 

If x = (q,p) G R 2 , and x' = (q',p') G R 2 , we write x A x' = qp' — pq'. Let A(t) = 
(a(t) ,b(t)) G R 2 , t G R, be any path, and as a special case of Proposition 23, suppose 
that G t is a linear function of x = (q,p): 

G t (q,p) = a(t)p-b(t)q = A(t) Ax (95) 

We get: 

H' (x, t) = ^A^x- ^A^j +H(x-A(t),t) (96) 

Proof. We have (^§f) = §q = ^ Ax, and since G t is linear, X G = X dasS)G = {G t , ■}■ 

Then X Gt (f*) = {G t , #} = f & - af t = § A A, so R = § A (x - \A). Also we have 
exp (— X Gt ) H = exp (— X c i asS) G t ) H = H (q — a(t) ,p — b (t)), i.e. the translated function, 
as transformed by the classical canonical transformation. □ 



D.2.4 Symplectic linear transformations 

Suppose that G t (q,p) is quadratic in q,p variables, with time-dependant coefficients. Then 
except from the quadratic terms who may change, the transformed Hamiltonian H' given 
in Proposition 23 is just the classical transform of H: 

H' = exp (—X c i asS) Ct) H + quadratic terms (97) 

Proof. Again Xc t = X c i aS s,Gt = {Gt,-}- Also, for any quadratic homogeneous functions 
A, B, then {A, B} is also quadratic. We deduce that X k lassGt (^) and R are quadratic 
homogeneous function of (q,p). Then Eq.(92) gives Eq.(97). □ 
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D.2.5 Time-dependant transformations of a special kind 



When dealing with next order terms, we will only need symbols G t of the simple factorized 
form: 

G t (x)=g(t)G(x) 

Then 

'dGt 



H' 



Of 



exp(-X Gt )H 



(98) 



since X Gt (fi) = g (t) g' (t) X G (d) = 0. 
D.3 Pre-Normalization 

Consider a given point Xo = (qo,po) G M 2 at time t — 0, and x (t) = (q (t) ,p (t)) = M jXo 
the trajectory passing through it. See figure 10. We will sometimes use the condensed 
notation x = {x (t) , t) e M 3 . 




Figure 10: Trajectory x(t). 



During Pre-normalization, we will perform successively two canonical transformations 
(a translation and a linear symplectic transformation) in order to simplify the Taylor 
expansion of H (q,p, t) along the trajectory x (t), up to degree two in (q,p). Higher degrees 
will be treated in a second stage called the normalization process. We will use the notation 
x' = (q',p') for the new coordinates and (q',p') for the transformed Hamiltonian, which 
depends on time through x = {x (t) , t). 

Notice that if F (x) — F (x, t) is a periodic function with respect to t, then it fulfills 
the obvious relation: 

F(x(t),t) t=2 = F((Mx)(t),t) t=0 (99) 
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D.3.1 Translation 

We decide to map the trajectory x (t) = (q (t) ,p (t)) onto the origin x 1 = (q',p') = 0, using 
a time-dependant translation, generated by Eq.(95), with A (t) = —x(t). From Eq.(96), 
the value of the transformed Hamiltonian H' at the origin is then 

h' x ( x ' = o) = - a 

with 

1 dx < i 1 / dq dp\ . . , 

A ^2T t hx - H[x) = 2{ p i- q i)- H ^- (100) 

We recognize the time derivative of the classical action along the trajectory x (t) (see [25] 
for a geometric interpretation). From Eq.(96) and using classical Hamilton equations for 
x(t), one obtains that the derivative of H' x at the origin is zero: (d x >H') (x 1 = 0) = 0, as 
expected because the origin is now a fixed point. Higher derivatives are equal to derivatives 
of if at x{t): 



D.3.2 Linear Symplectic transformation 

Similarly to Eq.(12), let Q x G SL(2,M) be the symplectic matrix which transforms the 
canonical basis of M 2 to the basis (u x , s x ) at point x (t) and time t defined above. We use this 
canonical transformation in order to simplify the quadratic terms of the Hamiltonian. This 
transformation is generated by a quadratic function G, and from Eq.(97), the transformed 
Hamiltonian H' differs from H only by different quadratic terms. From Eq.(7) and Eq.(87), 
the new quadratic terms are 



D.3.3 Result of the pre-normalization process 

Up to now, we have chosen an initial point Xo € M 2 , and have performed successively a 
translation and a linear symplectic transform: 

so that the trajectory x (t) is mapped onto the origin, and the new Hamiltonian (q',p') 
has a Taylor expansion in variables x' = (q',p') at the origin of the form: 



(#i) de9 >3(^0 = (#°W) (<?V) 

Terms of degree lower or equal to two are now in normal form, which means that the stable 
and unstable manifold are tangential to the principal axis, see figure 11. 
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t t 




Figure 11: The Pre-normalization transformation T x translates the trajectory to the 
origin, and transforms the unstable/stable directions onto the principal axis (q',p'). 

Eq.(91) with t = 0, t' = 2, gives 

M = %re,Mx M' XQ 'T pr l XQ 

where T pre ^ xo = T pr e,{x ,t=Q)- This relation gives Theorem 45 page 25 at the linear order 
J = 2 (we will explain below how to transform the time dependent term rj x q'p' into the 
time independent Ao,i )a . 9V, when treating resonant terms). 

Notice that the Hamiltonian (q ; , p') in pre-normal form we have obtained, is contin- 
uous with respect to x = (x,t), and periodic in t, as explained in Eq.(99). The action A x 
is not periodic in x with respect to the lattice r t , but the terms of higher degrees are. 

D.4 Normalization 

We go on with the normalization process, starting now from the Hamiltonian H x (q',p') 
given by the pre-normalization above, and which has a Taylor series in x' = (q 1 , p') at the 
origin of the form: 

# x (<?', p') = -a x + vWp' + VMxrtv'y' 13 (ioi) 

Z>0 a+(3>0 

In this expansion, every term is a periodic and continuous function of x = (x,t) G M 3 , 
(except for A x which is not periodic in x), and w a ,/3,/,x = 0, if / = and a + (3 < 2. We 
will transform Eq.(lOl) into a normal form, by successive iterations. At the first stage, we 
have no "quantum terms", i.e. Wa,/3,z, x = if / > 1. But during the normalization process 
such terms will possibly appear. We will check that due to the special form of the Moyal 
formula Eq.(90), only terms with even value of / will appear. 
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Notations: To simplify the notations, we will write rj (t) = 77( x (t),t), v a,p,i (t) = v a ,/3,i,(x(t),t) 
etc.. in order to show the time dependence along the trajectory x(t), and keep in mind 
that they also depend continuously and periodically on x and t. Let 

n d =a + (3 + 2leN (102) 

be the degree of the monomial h l q la p 113 . We will use the notation O (n) == O ((H, q' , p') n ) 
which means higher degree terms with respect to this graduation. 



D.4.1 Homological equation 

Suppose that we want to suppress the monomial term v aj pj (t) h!q' a p /l3 with degree n = 
a + P + 2l, from the Hamiltonian Eq.(lOl). For that purpose, we consider a time dependant 
canonical transformation generated by the symbol G x (t,q',p') = g x (t) h l q' a p ,l3 x WiP'), 
where g x (t) is a still unknown function, but supposed to be continuous and periodic with 
x = (q,p,t). The function xil'iP') i s a Co° cut-off function 7 with compact support such 
that x{q',P') = 1 for (q',p') E D 1 = {(q',p') / \q'\ 2 + \p'\ 2 < l}. From Eq.(98), the new 
Hamiltonian will be 

H'=(^)+ exp {-Xq) H 



\9t ) 



Consider (q',p') E D 1 . We have 



X G (H)(q',p') = X G (r)(t)q'p')+0(n + l) 
= rj (t) {G, q'p'} + 0(n + l) 
= rj (t) {p-a)g (t) tiq'y 3 + 0(n+l) 

So H' contains the monomial 

K {q\ p') = (vaju (t) + ^ + V (t) {pL-p)g {t)\ K l q^ 

and differs from H by this term and terms of higher orders only. We wish to suppress this 
term, by solving K = 0. This gives 

g x (t) = (c (X ) - f V x (s) e {a ~ P) So ^s' ds \ e -(«-/S) /* ^ds 



where C (x ) is a constant which depends only on x = M_ t (x) E M. 2 (i.e. the point of the 
trajectory at t — 0). The continuity in x is obvious, but we have to impose the periodicity 
relation Eq.(99), which writes: 

0*o(2) (2) = g(Mx ) (0) 



7 The cut-off function x ls necessary otherwise the canonical transformation generated by G is not 
defined other all (q',p') G M 2 . Think for example of G = Cp'q' 2 . 
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This gives a relation for the function C (xo), called the homological equation (this 
is usual in normal form calculations (see [3] chap 5, [45] p. 100): 



C (Mx ) = (C (x ) - P (xo)) e-^-^^o 



with 

^0 = / Vx(t)dt, 
JO 

which is the expanding rate along a piece of trajectory, and with 

r 2 

P(x ) = / v x (s)e( a -®foV*ds> ds 
Jo 

We have to check now that the above homological equation has a solution C (xo) con- 
tinuous with respect to io 6 I 2 . We have to consider two cases. 

The "non resonant case" a ^ (3: If a > [3, we can express C (x ) from the past image 
M~ 1 xq and iterate: 

C (x ) = (C (Af^o) - P (M-^o)) e~ ia - p)x ^o 

= (C (M- 2 x ) - P (M- 2 x )) e-M^^+^O _ P (M-^o) e^^" 1 ^ 

etc... 

Similarly, if a < (3, we use images of Xq in the future: 

C(x ) = P (x ) + C (Mx ) e -W- a)x *° 

= P (x ) + e~^- a)K o (P (Mx ) + C (M 2 x ) e -^~ a)XM ^) 

etc... The following Lemma shows that these iterations converge nicely. 

Lemma 12. (David DeLatte [19]). For a > (3, the solution of the homological equation 
Eg. (104) is 

oo 

C (x ) = ~J2 P ( M ~' x o) e - (Q - /3)A - t(xo) 

t=l 

which is a continuous function of xq, and with 

A_i(ar ) = X m*x , for t>\ 

-t<k<-l 

is the expanding rate along a piece of trajectory of length t. 
For a < (3, the solution 



(104) 
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+0O 

C (x ) = J^P (M'x ) e -0- a )M*o) 
t=o 

is continuous in Xq, with 

A 4 (xo) = X m*x , for t>\ 

0<k<t-l 

= 0, for t = 

Proof. From uniform hyperbolicity assumption, there exist \ m in > 0, t min > such that 
Vxo G T 2 , t > t m i n A t (xo) > \t\ \ m in- P (x) is a continuous function of x G T 2 , so bounded 
by \P {x)\ < P max . Now if ( x ) = - £^ (. . .) is the partial sum of C (x) as given 
above, with N G N, then Vx G T 2 , (C^ 1 ) (x) - (x)| < e~ x — N P max . This implies 
that (x) converges uniformly for iV — > oo, and that C (xo) is continuous in Xq. □ 

The "resonant case" a = (3: In that case we cannot solve the Homological equation, 
so we keep the monomial term of H . We can however make a small simplification of the 
Hamiltonian. Let us suppose that H has a resonant term rjij^h 1 (q'p'Y , where is time 
dependant because of x = (x (t) ,t). Let us define: 

dof f 2 

Az,j,0r o ) = / Vi,j,xdt 
Jo 

We decompose the resonant term into rjij^ = fji,j,* + \M,j,{x )- From Eq.(103), the first 
time-dependent term fjij^h 1 (q'p'Y can be eliminated after a time dependant canonical 
transformation given by g x (t) = — J* fjij^dt. There remains a time-independent term 
\\j,{x )ft (q'p'Y ■ In particular the resonant term (l,j) = (0, 0) is the action term Eq.(lOO), 
and gives 

A„,„, w = - [ A^dt = - £ (I (p§ -<!%)-H <*)) A (105) 

D.4.2 Result of the normalization process, and final proof of Theorem 7 

By successive iterations of the normalization process described above, we can discard every 
non resonant term of the Hamiltonian Eq.(57), up to a given degree. We have performed 
a sequence of time-dependant canonical transformations in a given order indexed by n = 
1,2,... n max , and x: 

% = T x (t)Q x U Gl x U G2x . . . U Gnmax ^ 

In this product, each term Ug„. x is generated by a bounded symbol G n>x (q' ,p' ,t) with 
compact support, equal to G n ^(t,q' ,p') = gi )a ,b,x.(t) h l q' a p'^ inside the disk Di of radius 
1. The resulting time-dependant Hamiltonian H^(q',p') is a total symbol which grows 
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quadratically at infinity and has a Taylor expansion in x' = (q',p f ) at the origin of the 
form: 

H*{q',p')= l - hu^WY +0(J + 1) 

2(l+j)<J 

i.e. the Taylor terms of degree less or equal to J are in normal form, with time- 
independent coefficients Xi,j,( x )- We can apply Eq.(91), with time t — 0, t' — 2, and deduce 
that 

where M JjXQ is generated by H xo{ $ = Op Wey i (H xo{t) ) on time two, i.e. M J)XQ = Mj, Xo , t =2, 
solution of dM J)XQ)t /dt = - (i/h) H Xo ^M J)XOtt , and Mj tXQtt=0 = Id. Also f x means T^ Xjt=0 ), 

and involves gi, a ,b,(x) = 9i,a,b,x ( £ = °)- 

In order to simplify the notations, we rescale time by a factor 1/2, so that t G [0, 1], 
and correspondingly we multiply Hx (t) by a factor two. We obtain finally: 

H x (q',p') = K x (q' lP ') + 0(J+l) (106) 

K x {q\p')= Yl \u^ l (q'py 

l+j<J/2 

The functions Xij t ( x ), and gi >a ,b,(x) are continuous with respect to x G M 2 . They are 
also periodic in x with respect to the lattice Z 2 , except for the action Ao,o,(x) ■ To deduce 
Eq.(25) of Theorem 7, we still have to compare M with the unitary map generated by the 
normal form part K x only, dropping out the terms of degree higher than J, and estimate 
the error. This is exactly what gives Theorem 20, page 49. Thus, we have finally proved 
Theorem 7. 



E Trace of a semi-classical hyperbolic normal form 

In this appendix, we give explicit semi-classical expressions for the trace and matrix ele- 
ments of a hyperbolic normal form. They are already obtained in [43], but we adapt them 
to the present context. Consider the Hilbert space H p i an e = L 2 (M), and the quadratic 
normal form Hamiltonian: 

K Q = Op Wey i (qp) = ^(pq + qp) = hi 

with 

t .1 / d d 



% 2 \dq^ ^ dq 



Let 

def 



k d ^ hflU ( 107 ) 

i,j>o 
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where jlij G E, but /io,i > 0. In this paper, the operator K appears in Eq.(60) and Eq.(61), 
as a result of the (post) normalization procedure along a classical trajectory. Remind that 
fil,j = if I < 1- The operator K generates a unitary propagator operator after time t: 

N (t) = exp (-itK/h\ (108) 
The truncation operators P a has been defined in Eq.(21). 



Proposition 24. For any < a < 1/2, one has a semi-classical expression: 

Tr (P a N (i) P a ) = T sem i + O (109) 

T semi d ^ f exp (-i^) exp (-it/xi,„) ^~ 
a semi- classical series 

s = i + J2 hSE s = l + hE x + h 2 E 2 + ... 

S>1 

where E s depends on t and jlij , with (I + j) < s + l, and are given explicitly in the proof. 
We have the control: 

\E S \ < t s E maXtS (HO) 

where E maxs does not depend on t. 



We prove this Proposition below. The idea and techniques of the proof are taken from 
[43] and related work cited in this paper. In this Proposition, the time t is fixed with 
respect to h~, but the result remains valid for \t\ = O (log (1/h)), because the semiclassical 
expansion remains unchanged. 

E.l A first useful formal asymptotic formula 

dcf A { 

Define /i (X) — Y2i ,->o 'X 3 ■ From Eq. (107) and Eq. (108), one can write N t = exp ( —j^t/j, 
Define: 

r (X) d ^ f y, (X) - (^,0 + h^o + lM>,iX) = ^ fa^X* 

l+j>2 

This decomposition is natural because for a given /6l, the expression: 

S = ~r (hi) = -it hj hl+j ~ lp = ° ( h ) (HI) 

l+j>2 
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is an asymptotic series in h: 

S = J2 hS S s (112) 

S>1 

with 

S s = (^-it^2fl(i =s -. j+ i)jP j , s > 1 

The first terms are explicitely: 

Si = -it (/2 2 ,o + £1,1 J + h^I 2 ) 
S 2 = -it (/2 3)0 + /22,i/ + Ai,2^" 2 + /io,3^ 3 ) , 

Then 

exp (--^tfA {fH) \ = exp f-^*// ,oJ exp (-it/i lfi I) exp (-t/i ,i) £ {I) (113) 
with the series in h and 7: 

£ (J) = exp (S) = l + £ ^ ( J ) ( 114 ) 

S>1 

where i£ a (J) is a polynomial of degree 2s in J: 

E S (I) ^ /'..,/' 
3=0 

The first two terms i? s (J) are: 

E l (I) = S 1 = -it (/2 2 , + + ^o, 2 / 2 ) (115) 
E 2 (J) = S 2 + ^S? = - (itfao + ^ 2 /l| J - 1 + ?fa,ofa,i) (H6) 

- J 2 ^t/ii j2 + ^* 2 At 2 ,i + t 2 ji2fijiQ^ - I 3 (it/2o,3 + t 2 h,ifa,2) (117) 

- / 4 (^,2) • (118) 
etc.. (119) 

However, in this paper, due to the special choice of Weyl quantization, we have fiij = 
if I is odd. This simplifies: 

E 1 (I) = -it (£ 2 ,o + /io, 2 / 2 ) 

£ 2 (/) = - Q^2,oJ " I {^,1) ~ I 2 (^2,0/i0, 2 ) " / 3 (^0,3) " I A Q^,2j • 
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We get: 



with 



and 



N (t) = exp ( -H^y J exp (-i*/x 1>0 ) exp I - % -r ( hi ) 



iVo (t) == exp ( —it^j^-K ) = exp ( — ifil 



def _ , 



We have for any k G N, 



so for any formal series £ (I) in J, this gives: 



(120) 



£[i^)N (t) = £[I)N (t) 



It means that we may substitute / <=?■ [i-jj- ) in the formal series. We deduce the formal 



result 



N (t) = exp ( -it^f-J exp {-itfo) £ (i^ ) N (t) 



:i2D 



Remark: In the right hand side of Eq.(121), only N (t) is an operator. 



E.2 Semi-classical expression of the Trace 

We first compute the Trace of Eq.(120) in the linear case. From Eq.(36), the functional 
T [ Nq) is we ll defined. It's value is 



T/=T(iV (t) 



2sinh(f) 



a* — •Vi^ 



(122) 



Remark: it can also be written T M = 1/ y |det (N (t) — I)\ where N (t) e SL (2, R) is the 
classical map. 

Proof. From its Schwartz kernel (x'\N (t) \x) = e^^ 2 5 (x' — e^x), and a suitable regulariza- 



tion, one has T (n (f)) = e» /2 J dx5 (x - e^x) = e^ 2 -^- 



2sinh(f 



□ 
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In Lemma 1 page 23, we showed that the value T M comes micro-locally from the origin: 

Tr (P a N (t)P a ) =T M + 0(/H 



From Eq.(114), £ T M = 1 + £ aJ h s E sJ T M = T,£ with 

2s 

s = i+ hSE s, E s = J2 E ^ lW 

s j=0 

and 

jCj) d ^ f ^ rfiT M 

In other words, we have substituted P by 7^ in the series Eq.(114). Finally, we deduce 
from Eq.(121) that: 

Tr (p a iV (i) P Q ) = exp (-H^f-J exp (-z^i, ) £ T M + O (ft 00 ) 

= exp ("i*^J exp (-itjl 1>0 ) £T^ + (ft 00 ) 

We have obtained Eq.(109). We can calculate the first terms 7^. With s = f sinh (|) 
and c = f cosh (^) , we obtain: 

J« = -^, 1(^-^(2^-^) (123) 

7® = - A (-6c 3 + 5cs 2 ) , 7< 4 ) = J-, (24c 4 - 28s 2 c 2 + 5s 4 ) . 
8s 3 V ; 16s 4 V ; 

Expressions for large values of t: If t S> 1, then /i = A ,it 3> 1, and 

= e-/ 2 (1 - e-)- 1 = 5>p (-n (k + 0) 



thus 



7« =(,r'4 ^=H) j ^E( A; +0 Jex p(-^( A; +0) 

w+x:fHy-H)V 



2 
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and we deduce in particular that 1^ is uniformly bounded for t E [1, +oo[. From Eq.(114), 
we deduce that 

\E I < t s E 

\ J -'s\ — - LJ max,s 

where E maX)S does not depend on t. 
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